ApJ, Accepted June 28, 1999 

Preprint typeset using lAT^^X style emulateapj 



DIFFERENTIAL GALAXY EVOLUTION IN CLUSTER AND FIELD GALAXIES AT Z « 0.3 

Michael L. Balogh^ Simon L. Morris^'^, 

H. K. C. Yee^''^, R.G. Carlberg^'^, and Erica Ellingson*''^ 

ApJ, Accepted June 28, 1999 

ABSTRACT 

We measure spectral indices for 1823 galaxies in the CNOCl sample of fifteen X-ray luminous clusters 

at 0.18 < z < 0.55, to investigate the mechanisms responsible for differential evolution between the galaxy 

Q> \ cluster and field environments. The radial trends of D4000, Wo{H5) and Wo{OII) are all consistent with 

0^ ■ an age sequence, in the sense that the last episode of star formation occurred more recently in galaxies 

0^ ' farthest from the cluster center. Throughout the cluster environment, galaxies show evidence for older 

stellar populations than field galaxies; they have weaker Wo{OII) and Wo{H5) lines, and stronger D4000 

5;^ ■ indices. 

^ ' From our primary sample of 1413 galaxies, statistically corrected for incompleteness and selection 

effects, we identify a sample of K+A galaxies, which have strong 115 absorption lines [Wo {H6) > 5 A) 
0^ ' but no [O II] emission {Wo{OII) < 5 A), perhaps indicative of recently terminated star formation. The 

CN \ observed fraction of 4.4 ± 0.7% in the cluster sample is an overestimate due to a systematic effect which 

results from the large uncertainties on individual spectral index measurements. Corrected for this bias, 
we estimate that K+A galaxies make up only 2.1±0.7% of the cluster sample, and 0.1 ± 0.7% of the 
field. From the subsample of galaxies more luminous than Mr — —18.8 + 51og/i, which is statistically 
representative of a complete sample to this limit, the corrected fraction of K+A galaxie s is 1.5+0.8% in 



> 

o 

—1. , the c luster, and 1.2+0.8% in the field. Compared with the z « 0.1 fraction of 0.30% (Zabludoff et al 

\^ 1996|), the fraction of K+A galaxies in the CNOCl field sample is greater by perhaps a factor of four. 



^ 



^ 



f~^ ' but with only la significance; no further evolution of this fraction is detectable over our redshift range. 

QV i We compare our data with the results of PECASE and GISSEL96 spectrophotometric models and 

0> ■ conclude, from the relative fractions of red and blue galaxies with no [OIIJA3727 emission and strong 

Hi5 absorption, that up to 1.9 + 0.8% of the cluster population may have had its star formation recently 

Q^. truncated without a starburst. However, this is still not significantly greater than the fraction of such 

I ' galaxies in the field, 3.1+1.0%. Furthermore, we do not detect an excess of cluster galaxies that have 

O \ unambiguously undergone a starburst within the last 1 Gyr. In fact, at 6.3+2.1%, the A+em galaxies 

that Poggianti et al. (1999) have recently suggested are dusty starbursts are twice as common in the 

^ ' field as in the cluster environment. 

Our results imply that these cluster environments are not responsible for inducing starbursts; thus, the 
increase in cluster blue galaxy fraction with redshift may not be a strictly cluster-specific phenomenon. 
We suggest that the truncation of star formation in clusters may largely be a gradual process, perhaps 
due to the exhaustion of gas in the galactic disk over fairly long timescales; in this case differential 
C^ ' evolution may result because field galaxies can refuel their disk with gas from an extended halo, thus 

regenerating star formation, while cluster galaxies may not have such a halo and so continue to evolve 
passively. 
Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: stellar content 
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is dominated in numbers by dwarf spheroidal galaxies and 
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in ma is by early type galaxies (E, SO, early spirals). There 1982 
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is strc ng theoretical (e.,e,., |Gott fc Gunn 1972|, [Antonuccio 
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Colafrancesco 1994| , Van Haarlem fc van de Wey 
1993 , Ibrman 1998p and observationa l (e.g., Zablud 



Franx 1993, Henriksen & Jones 1996) evidence that 



much of the cluster galaxy population has been built up by 
the accretion of galaxies from the surrounding, low den- 
sity, field environment, which implies that the rich clus- 
ter er Lvironmcnt has affected a morphological change in 



Bothun fc Dressier 1986 , Balsara, Livio fc O'Dea 
1994). If only gas in an extended, diffuse halo is stripped, 
star formation may be allowed to continue by consuming 
the remaining disk gas but, without infall to replenish this 
suppl y, star formation will die out on t imescales of a few 
Gyr ( p^jarson, Tinsley fc Caldwell 1980 ). Interaction with 



the ICM near the virial radius, where the ICM density 

is not high enough to completely strip the gala xy of its 

gas, may indu c e intense bur s ts of star formation (iDresslei 

these galaxies. A clear correlation between galaxy mor- fc Gunn 1983 , Evrard 1991, Gavazzi fc Jaffe 19871) which 



phology with local density ( pressler 198C ^ Dressier et al 



1997^ and cluster-centric radTus (IWhitmorc, Gilmore 



Jones 1993) has been observed. The well-known lack of 



emis sion line galaxies in clusters QOsterbrock 1960|7 Gisleij 



use up the remaining gas supply. Such bursts may also 
be induc ed by interactions with nearby galaxies (Moore et 
al. 1996| ,1998, Fujita 199^ ) or, near the clust er center, by 
interactions with the tidal field of the cluster (|B^rd fc Val' 



showr to correlate with cluster-centric radius (Balogh et 
al. 1997) and to be not e ntirely accounted for by the " 
morp holog y-radius relation (Moss fc Whittle 1993, Balogh] 
et al. 1998). Recently, Abraham et al. (1996) and Morris 
et al. (1998) showed that these various effects could be in- 
terpreted as an "age" gradient, in the sense that galaxies 
farther from the cluster center formed part of their stellar 
component more recently than the central galaxies. 

Another key piece of evidence that cluster galaxies are 
evolving strongly is the observation that the fraction of 
blue, spiral galaxies in cluster s incre ases with redshift 
(e.g. , [Butcher fc Oemler 197 j, 1984, [Dressier fc Gunn 



1998) 



19781 Dressier, Thompson fc Shectman 1985|) has also been tonen 1990| , |Valluri 1993| , [Henriksen fc Byrd 1996| , [Fujita 



1983, pressler et al. 1994, Couch et al. 1994, Rakos fc 



Sch onlibert 1995|, JYee ctTlJ. 1995[ , Pouch et al." 1998[ , but 
see pmail et al. 1998 ). This is usually referred to as the 
Butcher-Oemler (B-0) effect. To avoid possible misun- 
derstanding, we will always refer to the B-0 effect as this 
exclusively, a separate effect from the above problem of dif- 
ferential evolution between the cluster and field, though 
the two phenomena may be related. It has been shown 
that the star formation rate in the field also increases 
with 
al. 1 



edshift (iBroadhurst, Elhs fc Shanks 1988[, [Lill~ 



)95|, IRowan-Robinson et al. 1997|, [Tresse fc Mad 



5rp 



Understanding the mechanisms responsible for this 
change in star formation history is necessary to appropri- 
ately correct the CNO Cl universal mass density estimate 
(Carlberg et al. 1996) for the difference in mass-to-light 
(M/L) ratio in clusters, compared with the field. The cor- 
rection applied in Carlberg, Yee fc EUingson (1997) could 
be compromised if a large number of cluster galaxies un- 
derwent strong starbursts which do not occur in the field 
population. 

Dressier fc Gunn (1983) identified a galaxy population 
which they later termed E+A galaxies, that have strong 
Balmer lines but no detectable emission lines. These galax- 
ies (which we will hereafter term K-l-A to be consistent 
with current nomenclature) appear to have recently ter- 
minated star formation, within the ^1 Gyr prior to ob- 
servation. As such, they play an important part in under- 
standing the evolutionary history of galaxy populations. 
Evidence for these (and related) galaxies have been found 
both locally ( [CaldweU et al. 1993[, [Zabludoff et al. 1996[) 
and at moderate redshift (e. g., [ Sharpies et al. 1985 , Lav- 



ery fc Henry 198q, pouch fc Shar 



lis fc Couch 19881, OBroadhurst, E 



pies 1987, MacLarcn, El 



lis fc Shanks 1988, Fab 



dox 1 J98I, [Lin et al. 19991); thus, the B O effect as we 



ricant et al. 1991|, [Abraham et al. 1996 , Fisher et al 



consic cr it here may not be an exclusively cluster-specific 
effect and cannot naively be used as an explanation of the 
different populations which inhabit cluster and field en- 
vironments. Furthermore, it is not clear that the B-0 
effect can be directly interpreted as an evolutionary effect; 



1998[ , [Morris et al. 1998[ , pressler et al. 1999[ ), in both 
cluster and field environments. Many authors have in- 
terpreted the presence of these galaxies as evidence that 
a large percentage of cl uster galaxies recently u nderwent 
strong starbur s ts (e.g.. Dressier fc Gunn 1983 ^ Couch 



Kauff: nann (1995) has shown that rich clusters at z sa 0.4 Sharpies 1987 j_ Barger et al. 1996 ^ Poggianti fc Barbara 



1996, [Couch etal. 1998[, poggianti et al. 1999[). However, 



are not the direct progenitors of low redshift, compara- 
bly rich clusters. More recently, Andreon fc Ettori (1999) 
have shown that the higher redshift Butcher-Oemler clus- 
ters have higher X-ray luminosities than their low redshift 
counterparts. Since fair samples of clusters generally show 
either no evolution, or even n egative evolution in the X- 
ray l uminosity function (e.g., Henry et al. 1992|, [Collin j 



et al. [ 1997[, (vikhhnin et al. 1998URosati et al. 1998]), 



this confirms that the high redshift Butcher-Oemler clus- 
ters may not be the evolutionary predecessors of the low 
redshift clusters. 

Various mechanisms may be responsible for the trans- 
formation of a star-forming field galaxy to a passively 
evolving cluster member. An infalling galaxy which passes 
through or near the center of the cluster may have much of 
its gas content stripped by the hot intra-cluster medium 
(ICM), following which star formation will cease (Nulsen 



it is still debatable whether or not starburst and post- 
starburst galaxies observed in clusters are unique to these 
environments. Furthermore, other authors have claimed 
that most K-fA galaxies can largely be explained as the 
result of truncated star formation, without a large, initial 
burst (Newberry, Boroson fc Kirshner 1990| , Abraham et 
al. 199(: , Morris et al. 1998). It is necessary both to distin- 



guish between these two scenarios, and to compare their 
relative importance in the cluster and field environments, 
to determine which physical processes are plausible expla- 
nations for the differential evolution of cluster and field 
galaxies, and to determine how these processes will affect 
the M/L ratio in these environments. 

The CNOCl survey (Yee, EUingson fc Carlberg 1996) 
provides an excellent sample of galaxy spectra for which a 
reliable estimate of the star formation history for clusters 
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at z w 0.3 can be determined. This is a result of (1) a large 
sample size (over 2000 spectra); (2) well understood, and 
properly correctable, selection effects; (3) an identically 
selected field galaxy sample over the same redshift range; 

(4) objective measurements of spectral indices with well 
understood and empirically calibrated uncertainties; and 

(5) consistent comparison between data and spectropho- 
tometric models. From these data, we will attempt to 
specifically address the following questions. First, what 
fraction of galaxies in the CNOCl sample, regardless of 
environment, had significant star formation truncated in 
the last billion years or so? If this is truly larger than the 
low rddshift fraction, it suggests a link with the B-0 ef- 
fect. Secondly, how does the abundance ot these galaxies 
in the field compare with the abundance in clusters? An 
overabundance in clusters would suggest that the mech- 
anism responsible for generating these galaxies may be 
what leads to the differential evolution of cluster and field 
galaxies. Finally, how is the star formation in these galax- 
ies truncated? We will use spectrophotometric models to 
try to distinguish between a) truncation following a short 
burst; or, b) abrupt truncation without a burst. 

The paper is organized as follows. In §0 we define our 
data sample and spectral index measurements. In particu- 
lar, we discuss the details of our sample selection, statisti- 
cal corrections, magnitude limits, and index uncertainties. 
Interpretation of our lin e indices is made based on a com- 
parison with PEGASE ( Fioc fc Rocca-Volmerange 1997 ) 
spectrophotometric models, described in § |2.4 Our galaxy 
classification scheme is presented in |3, and is based on 
the line indices of local galaxies and results of the PE- 
GASE models (detailed in the Appendix) . We present our 
results in various forms in §0, and take care to explore 
systematic effects which result from the large uncertain- 
ties on our line indices. In §o we compare our results with 
the work of Zabludoff et al. (1996), Barger et al. (1996), 
Dressier et al. (1999) and others, and show that there is 
some consistency of results, though considerable variety 
in interpretation. We discuss some of the implications of 
these results in §pl and summarize our findings in m 

2. OBSERVATIONS AND MEASUREMENTS 
2.1. A Review of the Full CNOCl Sample 

The CNOCl cluster sample consists of CFHT MOS im- 
ages (Gunn g and r) and spectra of galaxies in the fields 
of fifteen^ X-ray luminous clusters in the redshift range 
0.18 < z < 0.55 selected from the EMSS survey"^ ( pioia 
& Luppino 1994). The observational strategy and full de- 
tails of the survey can be found in Yee et al. (1996). In 
particular, Yee et al. carefully consider and detail the se- 
lection effects of this sample, and describe how these are 
corrected for. For convenience, we repeat here some of the 
details in that paper relevant to the present analysis. 

The spectra are obtained with the O300 grism in place, 
which results in a dispersion of 3.45A per pixel and cov- 

^Omiliing cluster E0906+11, for which a velocity dispersion could not be computed due to an apparent double component structure(Carlberg 
et al. :996). 



ers the rest-frame wavelength range from (approximately) 
[OII]A3727 to the G-band (A « 4300A) in most galaxies. 
Slits are 1'.'5 wide, which results in a spectral resolution 
of about 16. 5A FWHM. Images and spectra of galaxies 
in each cluster are obtained from either one, three or five 
MOS fields mosaiced east-west or north-south, which pro- 
vide non-uniform coverage as far out as 1-2 R200 in pro- 
jected distance for most clusters, where i?2oo is the radius 
(« 1.2 h"^ Mpc) at which the mean interior mass den- 
sity is equal to 200 times the critical density, and within 
which it is expected that the galaxies are in virial equi- 
libriu m ( Gott fc Gunn 1972 ; Crone, Evrard fc Richstonc 
1994). For flo — 0.2, the cluster virial radius is equal to ap- 
proximately 1.2_R200- The final catalog contains redshifts 
for about 2500 galaxies, including field galaxies located in 
front of and behind each cluster. 

In order to average over effects of non-sphericity, we 
will combine all fifteen clusters to produce a single sam- 
ple, hereafter referred to as the cluster. Cluster member- 
ship is determined based on the observed radial velocity 
difference from the brightest cluster galaxy (BCG*). We 
use the mass model of Carlberg et al. (1997) to determine 
the (projected) radial dependence of the cluster velocity 
dispersion, cT{r), from the average, measured dispersion, 
(Ti. In Figure |l] we plot the galaxy velocities, normalized 
to (Ti, against their projected radius from the BCG, nor- 
malized to -R200- The 3cr and 6ct contours of the mass 
model are shown as the solid and dashed lines, respec- 
tively. Galaxies with normalized velocities less than 3cr(r) 
are considered cluster members (squares) , while those with 
normalized velocities greater than 6cr(r) comprise our field 
sample (circles). The population with intermediate veloc- 
ities is represented by the solid triangles, and are consid- 
ered "near-field" . The extent of the vertical axis in this 
figure is limited for clarity; hence only a small fraction of 
the total field population is displayed. The cluster popu- 
lation appears well separated from the field in this figure, 
and exhibits a rising velocity dispersion toward the clus- 
ter center. We will only include galaxies within 2i?200 in 
our cluster sample, as larger radii are much more sparsely 
sampled. 

By extracting the spectra with standard IRAF^ rou- 
tines, an error vector is created which includes Poisson 
errors, read noise and sky subtraction uncertainties. This 
allows a good determination of both the signal-to-noise 
(S/N) ratio of the data, and also reliable errors on any 
measured quantities (see §2^). The S/N ratio is defined 
as the mean value of the flux per pixel divided by the er- 
ror per pixel, in the wavelength range 4050 < A/A< 4250. 
The resulting S/N distributions for the cluster and field 
samples are shown in Figure g. The mode of the distribu- 
tion is at S/N«7 per pixel (or 15 per resolution element), 
and 94% of spectra have S/N^ 3. 

The entire CNOCl data sample, including the raw data 
and all measured quantities included in this paper, will 



^Rycept Ahell ^.SQO . which is not in the EMSS sample, but of comparable richness and X-ray luminosity to the other CNOCl clusters 
(Abraham et al. 1996). 



"Except for cluster MS 0451.5-1-0250, for which no redshift is available for the BCG. The velocities and positions are measured relative to 
the mean for this cluster. 

^IRAF is distributed by the National Optical Astronomy Observatories which is operated by AURA Inc. under contract with NSF. 
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Fig. 1. — This figure shows how we define our cluster and field samples. The data plotted are radial velocity differences 
from the BCG, normalized to the cluster velocity dispersion, cti, against projected radii from the BCG, in units of i?200- 
The solid and dashed lines are the 3cr(r) and 6a{r) contours, respectively, of the cluster mass model of Carlberg et al. 
(1997). We define cluster members {squares) as those which lie within 3a{r), and field galaxies (circles) as those with 
velocities greater than 6a{r). The intermediate population (solid triangles) are termed "near-field". The extent of the 
vertical axis is limited to clearly show the cluster population; only a small fraction of the field sample is shown in this 
figure. 
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Fig. 2. — The distribution of signal-to-noise per pixel is shown for the cluster (solid line) and field (dashed line) samples. 
The number of galaxies in each bin is weighted as described in § |2.3.1 and, thus, represents the number of galaxies in the 
photometric sample, though the distribution is computed from the spectroscopic sample. 
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soon be available from the Canadi an Astronomy Data 
Centre (http://cadcwww.liia.nrc.ca/). All cosmologically 
dependent calculations in this paper are made assuming 
go = 0.1, A = 0. The Hubble constant is parameterized as 
usual by /i = i7o/100 km s^^ Mpc"^ 

2.2. Spectral Index Definitions 

We will consider three spectral indices, tabulated in Ta- 
ble H, which are automatically measured by an IRAF IM- 
PORT routine available from the first author on request. 
The break strength at 4000A (D4000) is defined as the 
ratio of the flux in the red continuum to that in the blue 
continuum. We have chosen a much narr ower definition o f 
this break than the standard definition ( Hamilton 19851 ). 
There are two principal reasons for this. First, the un- 
certainty in the standard index value, assuming Poisson 
noise, was found to be a poor representation of the differ- 
ence between multiple measurement s of t he same galaxy 
(this analysis is discussed further in §2.2.1). Reducing the 
width of the continuum definitions significantly improved 
this correspondence. Secondly, the new index is much less 
s ensiti ve to reddening effects, which may be important (see 



4.31D . 

There are considerable advantages in using D4000 in- 
stead of the available g-r color. Most importantly, com- 
puting a rest frame color from the observed photometry 
requires calculating K-corrections from the very models 
we will be using to determine star formation histories and, 
thus, introduces a circular argument. Furthermore, the 
galaxy spectra are obtained from slits which only sam- 
ple the central regions of the galaxy, while the photome- 
try represents the integrated light from the whole galaxy. 
Thus, the photometric colors and observed spectral fea- 
tures may not arise from the same environments. (How- 
ever, at z = 0.3, the slit width of 1'.'5 corresponds to a 
linear size of AAh^^ kpc, which is a sizable region). In 
addition, especially for our unusually narrow definition of 
D4000, this index is very reddening insensitive. The prin- 
cipal disadvantage is that the measurement uncertainties 
in the D4000 measurements are generally larger than those 
of the g-r colors. 

The rest frame equivalent widths of the H5 Balmer 
absorption line, Wo{Il5), and the [OIIJA3727 emission 
line, Wo{OII)^ were automatically computed by summing 
the observed flux (accounting for partial pixels) above 
(Wo(0//)) or below {Wo{H5)) the continuum level, which 
itself is estimated by fitting a straight line to the flux in 
the continuum regions. The continuum data are weighted 
by the inverse square of the noise (from the IRAF noise 
vector), and the fit is then constructed by weighted linear 
regression. The blue and red continuum regions are de- 
fined in Table Ol the line flux is summed over the region 
defined "Line" in the same table. Note that the Wo{OII) 
index is positive when the line is in emission, while the 
WoiHS) index is positive when the line is in absorption. 

It must be noted that the indices defined above (or, in- 
deed, any such indices) are not necessarily related simply 
or directly to physical quantities. The value of the Wo{II6) 
index is very sensitive to the line and continuum defini- 
tions, and also on the spectral resolution and sampling. 
In particular, negative values of Wo{H5) do not always 
indicate emission, since absorption lines in the continuum 



regions (especially in galaxies with old stellar populations) 
may result in negative Wo{IId) index values. Thus, cau- 
tion must be taken when comparing these measurements 
with those of other authors; in fact, the index definition 
here differs slightly from that used by our own collabora- 
tion in Abraham et al. (1996). 

2.2.1. Index Uncertainties 

Errors on the spectral indices are determined from the 
noise vector generated from the optimal IRAF extraction. 
The D400Q index is simply a flux ratio and, hence, the er- 
ror is determined from standard propagation techniques. 
Errors on equivalent width measurements are computed 
from equation A8 in Bohlin et al. (1983); this makes use of 
the uncertainties of the weighted, linear continuum fit, and 
the uncertainties of each pixel defined as comprising the 
line itself. Often, spectral regions contaminated by bright 
night sky line emission are very noisy, or interpolated over 
by hand. Since the computed error vector reaches very 
large values at the wavelengths of night sky lines, this ef- 
fect will be reflected in the index uncertainties. Excluding 
from the sample all galaxies in which either of the lines 
at 5577 A, 5890 A, or 6300 A lands within the index def- 
initions of HS or [O II] does not significantly change our 
main results. 

Since spectra were obtained with up to three masks 
for each cluster and, also, since adjacent fields overlap by 
about 15", some galaxies were deliberately observed more 
than once. This allows us to test the reliability of our 
error estimate, by comparing the difference between two 
measurements (xi — X2) with the quadrature sum of their 
errors, \/<j\ -\- a 2- In Figure y we plot the distribution of 

the ratio of these two numbers, e = (xi — x-i)! \/ a\ + a\, 
for each of the three indices. If the true errors are Gaus- 
sian distributed with a variance given by our error es- 
timate, the distribution of e should be Gaussian with a 
mean of zero and variance of unity; this is shown as the 
solid curve for reference. All of the e distributions have 
a variance greater than 1, and are inconsistent with the 
solid curves. This implies that our raw uncertainties are 
underestimated, which may reflect the small but inevitable 
systematic error resulting from the subtraction of incorrect 
sky levels in some cases. The 3tT clipped variance of each 
e distribution is 1.44 for WoiOII), 1.42 for Woiji^) and 
1.96 for D4000. The distributions are better represented 
by a Gaussian of the corresponding, wider variance, shown 
as the dashed lines. To compensate for this effect, we mul- 
tiplied all error estimates by the appropriate value. The 
final index value adopted for multiply observed galaxies 
is the average, weighted by the square of the error, of all 
independent measurements (up to three). 

2.3. The Data Sample 

In order to carry out our analysis, statistical weights 
must be computed for all galaxies, to correct for the in- 
completeness of the spect roscopic sampling; this proce- 
dure is reviewed in §2.3.1. W c then proceed to define a 
luminosity-limited sample, in § 2.3.2| , which will be consid- 
ered for those analyses which are sensitive to such a limit. 
For the remainder of the analysis, for which the primary 
purpose is to compare properties in the cluster and field 
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Table 1 




Index 


Blue continuum 

(A) 


Line 

(A) 


Red continuum 

(A) 


D4000 

[OII]A3727 

H5 


3850-3950 
3653-3713 
4030-4082 


3713-3741 
4082-4122 


4000-4100 
3741-3801 
4122-4170 



Note. — Definitions of fine indices measured on the data and 
spectrophotometric models. The D4000 index is the ratio of the 
flux in the red continuum to that in the blue continuum. The other 
three indices are equivalent widths of the line over the wavelength 
range listed in Column 3, where the continuum level is estimated 
from the flux blueward (Column 2) and redward (Column 4) of 
the line. 




(xj-x,)/(a2 + a^)V2 



Fig. 3. — For multiply observed galaxies in the sample, we plot the ratio of the difference between two independent 
measurements of each index and the quadrature sum of the estimated index errors. If the true errors are normally dis- 
tributed and well represented by the estimated errors, the distribution will be Gaussian with a mean of zero and variance 
of unity, shown by the solid line. The distributions are wider than this Gaussian, but consistent with one of variance 1.44 
{Wo{OII)), 1.42 {Wo{H5)) and 1.96 (D4000), shown as the dashed lines. 
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samples, we define a larger sample which is not luminosity 
lim ited, b ut in which the lowest quality data are excluded, 
in § 2.3.3 . Unless otherwise specified, this latter is the data 
sample that will be discussed in the remainder of the pa- 
per. 



2.3.1. Data Wcight3 

Since spectra are not obtained for all of the galaxies 
observed photometrically, we must correct for this incom- 
pleteness and any selection effects that may arise as a re- 
sult. These effects, and the calculation of compensating 
statistical weights, are considered in detail in Yee et al. 
(1996); we review them briefly here. The main selection 
criterion is apparent magnitude; a smaller fraction of faint 
galaxies are observed spectroscopically, relative to brighter 
galaxies. The magnitude weight Wm compensates for this 
effect. A second order, geometric weight Wxy is computed 
which compensates for such effects as the undersampling of 
denser regions, and vignetting near the corners of the chip. 
Finally, a color weight Wc is computed to account for the 
fact that bluer galaxies are more likely to show emission 
lines, facilitating redshift determination; however, this is 
a small effect, and we have checked that the inclusion of 
this weight does not significantly affect any of the results 
discussed in the present paper. Both Wxy and Wc are nor- 
malized so that their mean is 1.0 for the full sample. In 
all of our statistical analysis, each galaxy in our sample is 
weighted by VFgpoc = Wm x Wxy x Wc to statistically cor- 
rect for these selection effects^°. We have also computed a 
"ring" weight, which is meant to correct for the fact that, 
due to the geometry of the cluster mosaics, the outer re- 
gions of clusters are less well sampled (in area) than the 
central regions. This correction is only important in the 
presence of radial gradients, if a globally averaged clus- 
ter quantity is sought. We will not use this weight in the 
present analysis, though we have verified that its use does 
not change any of the results (recall that we restrict the 
cluster sample to those galaxies within R < 2_R20o)- 

2.3.2. Luminosity Limited Sample 

We will first define a sample limited only in luminos- 
ity. Absolute r magnitudes (M^) are calculated from the 
photometry, and k-corrections are made based on the g-r 
colors and the model spectral energy distributions of Cole- 
man, Wu & Weedman (1980), convolved with the filter 
response function, for four, non-evolving spectral types 
(E/SO, Sbc, Scd and Im). We also correct M^ for redshift 
evolution by assuming galaxies brighten by a factor (l-l-z) 
(Lilly et al. 1998; Abraham et al. 1998); thus, we correct 
all luminosities to the corresponding z = values by di- 
viding them by (1+z), though this does not significantly 
affect our results. We chose an absolute magnitude limit 
of Mr — — 18.8 -f 51og/i, since this excludes most of the 
galaxies with large magnitude weights, Wm > 5; this limi t 
corresponds to 1.8 mag below M* in r ( King fc Ellis 1985| ). 
The absolute magnitude distribution of the cluster sample 
(weighted by Wgpcc) is compared with that of the field in 
Figure 0. There are 1125 galaxies in this sample (omitting 
the BCGs and cluster members beyond 2i?2oo)j including 

^''AU uncertainties in this paper that are based on vN Poisson 
consideration. 



710 cluster members, 343 field galaxies and 69 near field 
galaxies. 

The absolute magnitude distributions of the low (z < 
0.33) and high {z > 0.33) redshift galaxies in this sample 
are shown in Figure H. As a result of the fainter appa rent 
magnitude limit in the high redshift cluster samples (Yee 
et al. 1996| ), and the luminosity evolution correction, the 



two distributions are comparable. 



2.3.3. The Maxim,al Sam,ple 

For the purposes of comparing cluster and field galaxy 
populations, it is desirable to use as much of the data as 
possible, excluding only that of the poorest quality, and en- 
suring that the luminosity distributions of the cluster and 
field samples remain similar. From the 1823 galaxies for 
which both [OIIJA3727 and US lie within the observed red- 
shift range (excluding the BCGs and cluster members be- 
yond 2i?20o)i we first select the 1572 which have Wm < 5, 
so that the lowest quality data do not dominate our re- 
sults; this is equivalent to imposing a magnitude limit 
which varies from cluster to cluster. It is further desir- 
able to remove those galaxies from the sample which have 
extraordinarily large uncertainties in the measured spec- 
tral indices, which often result from poor subtraction of 
bright night sky lines. The distribution of the error es- 
timates for each index, scaled as discussed at the end of 
the previous subsection, are shown in Figure |6|. Based 
on these distributions, we have chosen to exclude an ad- 
ditional 159 galaxies with errors greater than 15A, 5A, 
and 0.5 for Wo{OLI), Wo{Hd) and D4000, respectively. 
This selection introduces another second order correction 
to Wspec, as galaxies with index errors exceeding our limit 
tend to be fainter. We calculate, in bins of absolute lu- 
minosity, the fraction of galaxies, weighted by M^spoc, re- 
moved from the sample by this selection. The remaining 
galaxies are weighted by the inverse of this number, which 
we will call Wom multiplied by M^spcc- W^crr varies from 1 
for galaxies with Mr < —19.8 + 5 log h, to ^ 2 for galaxies 
with Mr « —17.1 -I- 5 log ft,. The maximal sample consists 
of 1413 galaxies: 924 cluster, 407 field and 82 near field. 
The absolute magnitude distribution of the cluster sam- 
ple (weighted by W^poc x Worr) is compared with that of 
the field in Figure 0. Although the two distributions are 
similar, the cluster sample is biased toward less luminous 
galaxies, and the weights do not sufficiently correct the 
spectroscopic numbers (to be representative of the photo- 
metric sample) below Mr = —18.8-1- 51og/i, due to the 
removal of galaxies with Wm > 5. 

2.4. The PEGASE Model Parameters 



As discussed in § p.2[ , the Wo (HS) index is sensitive to its 
definition and the manner in which it is measured. Thus, 
its interpretation depends on comparing it with models in 
which the index is measured in an identical manner. It is 
not possibl e to use model results published by other au- 
thors (e.g.. Couch & Sharpies 1987, Barbaro & Poggianti 



1997) to interpret data unless it is shown that the indices 
dehned for those models are comparable to those defined 

statistics are computed from the unweighted number of galaxies under 
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Fig. 4. — The evolution and k-corrected, absolute Gunn r magnitude distributions of the clu ster ( solid line) and field 
{dashed line), luminosity limited galaxy samples. Each bin is weighted by H4pec, discussed in § ^.3.1 ; thus, the height of 
each bin represents the number of galaxies in the photometric sample, as determined from the spectroscopic sample. 
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maximal sample. This selection criteria is compensated for by a statistical weight, Werr, discussed in §2.3.3 
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for the data. In light of this, we will reconstruct models in 
§0 which have been explored by other authors, but based 
on our index definitions. 

The evolution of Wo {H6) and D4000 with star forma- 
tion history for solar metallicity galaxies i s obtained from 
the P EGASE pop ulation synthesis code (Fioc & Rocca-| 



range 1997), based on their UV to NIR spectral 



The PEGASE models predict results similar to 



Char- 



Volmc 
library. 

those of the updated GISSEL96 models (Bruzual 
lot 19B3), but provide more freedom in constructing spec- 
tra with complex star formation histories; in particular, 
different initial mass functions (IMFs) are readily imple- 
mented. The drawback is that the models are solar metal- 
licity; to investigate metallicity effects, we will consider 
the GISSEL96 models. The standard IMF adopted here is 
that of Salpeter (1955); the effec t of ad opting a truncated 



IMF model will be discussed in §4.3.1 



Spectral indices are determined from the model spec- 
tra in the same manner as from the CNOCl data; that 
is, the same code, with the same line definitions, are run 
on the spectra output by PEGASE. We have not included 
Balmer-filling due to nebular emission; thus, the model 
Wo{H5) indices represent an upper limit on the equiva- 
lent width one expects to observe when massive stars are 
present. 

Care was taken to ensure that the line indices we obtain 
from the models correspond to the line indices measured 
on the data; in particular, the Wo{Hd) index is quite sen- 
sitive to both the resolution and the sampling of the spec- 
tra. The (rest frame) resolution of the model spectra is 
lOA, comparable to the resolution of our data (about 15A 
at z=0.3). However, the model spectra are more coarsely 
sampled than the CNOCl data (every lOA instead of every 
2.6A, rest frame), so we resample them to 2.6A, interpo- 
lating using a fifth order polynomial fitting routine. We 
verified that each line index we measure on the CNOCl 
spectra is unchanged if we first smooth the spectra by 
averaging every four pixels, and then resample using our 
interpolation algorithm. 

3. GALAXY CLASSIFICATIONS 

There is, as yet, no clear consensus in the literature 
about what type of galaxy is called, for example, a "star- 
burst", "H(5-strong" , "post-starburst" , "E-fA/K-hA", 
etc. This is because the diagnostic indices are not all de- 
fined in the same way, the spectra used vary in terms of 
S/N and resolution, and the physical phenomena one is 
trying to isolate varies from author to author. In any case, 
such classifications are somewhat arbitrary, and there can 
be no objectively correct choice. For reasons discussed in 
§g, we have chosen not to use the most recent classifica- 
tion scheme (Dressier et al. 1999), and thus wish to state 
clearly how our classes are defined. These should be com- 
pared and contrasted not only with that paper, but also, 
for example, the differing definitions in Couch & Sharpies 
(1987), Barbaro & Poggianti (1996), Barger et al. (1996) 
and Morris et al. (1998). In most cases, the differences 
in classification are ones of detail, and refined precision of 
the evolving nomenclature. 

3.1. Determination of the HS Threshold 



One of the main goals of this paper is to identify a pop- 
ulation of galaxies which were forming stars a short time 
(< 1 Gyr) ago, but have since stopped. These must be 
distinguished from galaxies which are currently undergo- 
ing star formation, and those which have not done so for 
a long time. We will make this distinction partly based 
on the spectra of local galaxies (see § |3.2D , and partly on 
model predictions, which is the purpose of this subsection. 

In Figure S we show, as the dashed line, the evolution 
of the Wo{H6) index for a galaxy with a constant star 
formation rate. As long as massive stars are present, neb- 
ular emission filling of US will reduce the index by be - 
tween about lA and 3.5A (Barbaro & Poggianti 1997); 
thus, such a spectrum will not reach Wo{HS)>5A.. All nor- 
mal, star-forming galaxies are therefore expected to have 
Wo{HS)<^ 5 A, and nebular emission lines. This is con- 
firmed by the line indices of normal, spiral galaxies, as we 
show in § |3.2| . 

To reach Wo{HS)> 5A, the spectral light of a galaxy 
must be dominated by A- and early F-type stars. The 
presence of O- and B-type stars, which have weak intrin- 
sic as absorption and can dominate the galactic light, is 
the reason why the dashed line in Figure |8| does not reach 
Wo{HS)> 6.5A, even without considering nebular emis- 
sion contributions. For illustrative purposes, we show the 
evolution of Wo (HS) in Figure @ for a model galaxy with 
a constant star formation rate, but with no stars more 
massive than 2.7 Mq (corresponding to an AO star, Bohm- 
Vitense 1992). There will be no emission contribution to 
such a spectrum, which therefore has Wo{HS)> 5 A for up 
to 7 Gyr; after that the accumulation of less massive stars, 
which have weak Wo {HS) absorption, begins to contribute 
significantly to the light. Clearly the absence of OB-stars 
allows Wo{HS) to reach muc h higher levels. As we will 
discuss in more detail in §3.3, strong Wo{HS) indices may 
be observed following the abrupt termination of star for- 
mation, due to the almost immediate death of the O- and 
B-type stars. This can be visualized in Figure g as follows. 
A constantly star-forming galaxy will evolve at least ~ 1 
A below the dashed line; as soon as star-formation termi- 
nates, Wo{HS) will increase approximately to the height 
of the solid line (which is an A-star dominated spectrum) . 
For this reason, terminating star formation after 3 or 4 
Gyr results in only a weak increase in Wo {HS) , while ter- 
minating a short burst of a few hundred Myr results in a 
much larger increase. In either case, emission lines will be 
absent in the spectrum. 

Based on the results of these models, we adopt a thresh- 
old of Wo{HS)=5 A as an acceptable division between 
galaxies with normal rates of star formation, and those 
with recently truncated (or otherwise anomalous) star for- 
mation. This threshold is confirmed by the line indices of 
normal local galaxies, as discussed in the following subsec- 
tion, and may be lowered to Wo(HS)—3A for the reddest 
galaxies, as we demonstrate in § p.3| . 



3.2. Definitions Based on Wo{OII) and Wo{HS) 

We will classify the galaxies in our sample based prin- 
cipally on their position in the Wo{OII)-Wo{HS) plane 
relative to the PEGASE models, and to local galaxies. 
We show this plane in Figure ^. The data plotted are 



Balogh et al. 



11 




rime (Gyr] 



Fig. 8. — The evolution of Wo{HS) with time, as determined from PEGASE models. The dashed line traces the evolution 
of a galaxy with a constant star formation rate, neglecting nebular emission, which will reduce Wo{HS) by ~ 1 — 3.5 A. 
The solid line is the same model, but without stars more massive than 2.7Mq; there will be no emission contribution in 
this case. The horizontal line at Wo{HS)—5A is the lower Wo{HS) limit for our K+A galaxy definition. 



indices measured by running our code with our line defi- 
nitions on the local data of Kennicutt (1992b, open sym- 
bols) and Kinney et al. (1996, filled symbols). The dif- 
ferent symbols represent the different morphological clas- 
sifications as shown on the figure. The diamonds include 
any galaxy classified as peculiar, irregular or starburst; all 
other points correspond to "normal", local galaxies. The 
plane is divided into five regions, each of which we will 
discuss separately. 

• Passive Galaxy: This is a galaxy that is not 
currently undergoing significant star formation, and 



:hus evolves passively. Spectrally^ we identify such a 
galaxy as one with Wo{OII)< 5A, which is roughly 
our detection limit, and Wo{H6)< 5A. The local 
galaxies in this category are mostly elliptical, SO 
or early type spirals; the presence of these normal 
galaxies with Wo{Hd) as large as nearly 5A justifies 
the need for a threshold of at least this strength to 
clearly isolate galaxies with anomalous star forma- 
tion histories. 



• Star Forming (SF) Galaxy: This is a galaxy 
that has been undergoing significant star formation 
for at least several hundred million years. Unless 
heavily obscured by dust, the [OIIJA3727 emission 
line should be prominent, so Wo{OII)> 5 A. Further- 
more, models show that <Wo{H6)< 5A; the lower 
limit is determined from the results of Poggianti & 
Barbaro (1996) and Barbaro & Poggianti (1997), 
which show that, if star formation has been ongo- 
ing for at least several hundred million years, the 
stellar absorption will dominate any nebular emis- 
sion of this feature. Local galaxies that lie in this 



category are both late Spirals (Sc) and some galax- 
ies that have been classified "starburst", based on 
the presence of strong emission lines. 

Short Starburst (SSB): We define a short star- 
burst as a short lived (< 200 Myr), large increase 
in star formation rate. A galaxy that has had a 
high star formation rate for more that a few hun- 
dred Myr will be considered a SF and not a SSB 
galaxy. In this case, nebular emission will dominate 
the Wo{H5) feature and, hence, b oth 1^0(0//) and 



Wo{H5) will be emission features (Poggianti 
baro 199^, [Barbaro fc Poggianti 1997|). 



Bar- 



• K+A Galzixy: This type of galaxy is one 

with no detectable [OIIJA3727 emission (defined 
as Wo{OII)< 5A), but strong Hi5 absorption 
{Wo(H5)> 5A). Such a spectrum, which requires 
the presence of spectral type-A stars without more 
massive stars (which would give rise to [Oil] emis- 
sion lines), may be the result of a recently termi- 
nated episode of star formation, most likely follow- 
ing a short starburst (se e, for example. Dressier 
& Gunn (1983), and j |3.3| ). Alternatively, some of 
these galaxies may be extreme examples of the dust- 
obscured starburst scenario advocated by Poggianti 
et al. (1996) to describe the A-|-em galaxies (see be- 
low). This simple classification is based on only two 
line indices, and our K+A galaxies do not comprise a 
completely homogeneous class in terms of other spec- 
tral features, including the strength of other Balmer 
lines, or the Ca II H and K lines. This must be kept 
in mind when comparisons are made with other work 
in which these features are considered. 
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Fig. 9. — Local data of Kennicutt (1992b, open symbols) and Kinney et al. (1996, filled symbols) presented in the 
Wo{OII)-Wo{H5) plane. Positive values of Wo{OII) and negative values of Wo{HS) represent emission in these indices. 



The plane is divided into five regions according to the definitions in §3^: Star-forming (SF) galaxies are those which 
have undergone significant star formation for at least several hundred Myr, while star formation in short starburst (SSB) 
galaxies has begun within about 200 Myr. The origin of A+cm galaxies, which exhibit a strong A-star spectrum with 
emission lines, is uncertain; some may be dusty starbursts, or host an AGN. Of galaxies with Wo{OII)< 5 A, K+A galax- 
ies may have recently had their star formation truncated, while those galaxies with weaker Wo{HS) indices are evolving 
passively, with little or no star formation. 



• A+em Galaxy: Our final category consists of 
galaxies with Wo{OII)> 5A and Wo{H6)> SA. This 
corresponds to a galaxy with Balmer absorption lines 
as strong as those found only in A- and early F- 
type stars, but with [OIIJA3727 emission lines in- 
dicating the presence of more massive O- and B- 
type stars. There is not yet a clear interpretation of 
these spectra; the presence of OB-stars should keep 
Wo{H5)<5K due both to emission-filling and the 
weaker, intrinsic Wo{H6) absorption of these lumi- 
nous stars. It is possible that some of the observed 
emission lines are due to an active nucleus; the useful 
diagnostic emission lines (Ha,H/3,[N II], [O III]A5007) 
are generally redshifted out of our spectral range, so 
that we are unable to confirm this possibility. Alter- 
natively, Poggianti et al. (1999) suggest that these 
galaxies (which they call e(a)) may be dust -obscured 
starbursts; we discuss this point further in dq, and 
note that the same interpretation may hold for the 
previously discussed K-f A class of galaxy. The lone 
point in the Kennicutt (1992b) sample which lies in 
this region, shown on Figure |^, is NGC 3034, which 
is in fact a very dusty, star-for ming galaxy. It has 
been classified "starburst" (e.g., Rieke et al. 1993 ) 
based on its strong Ha emission, but has long been 
known to have an inte grated stellar sp ectrum anal- 
ogous to A-stars (e.g., Humason 1956); this galaxy 
would thus seem to be of the type advocated by Pog- 
gianti et al. (1999). 



3.3. Additional Definitions Based on D4-000 

To investigate the physical origin of the K-l-A galax- 
ies, we will compare our observations with the PEGASE 
spectrophotometric models to predict how Wo {H6) evolves 
with galaxy color, as determined by the D4000 index. This 
largely follows the method outlined by Couch & Sharpies 
(1987) and Barger et al. (1996), applied to the PEGASE 
models with our current index definitions. 

Barger et al. (1996) divided their sample into red and 
blue halves at B-R=2 (observed at z = 0.31). To convert 
this to an equivalent cut in D4000, we calculate observed- 
frame {z = 0.3) B-R and rest frame D4000 for two different 
star formation histories, using the PEGASE models. We 
show the results in Figure |l^; the models are described 
in the caption. From this figure, we adopt D4000=1.45 
as the value corresponding to B-R=2; this value is fairly 
independent (within 0.05) of the model considered. 

We now plot the output of five different models in the 
DAOOO-Wo{HS) plane, shown in Figure [III and described 
briefiy in the caption; we defer detailed discussion of these 
models to the Appendix. The local data plotted on this 
Figure are the same as those shown in Figure M. The 
plane is divided into five regions similar to those defined 
by Barger et al. (1996), with the help of the D4000=1.45 
division justified above. We have also renamed some of 
the classifications to correspond with our present termi- 
nology. The Passive, SF and SSB regions are expected to 
contain galaxies similar in star formation history to those 
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Fig. 10. — Observed frame B-R is plotted against rest frame D4000 for two different PEGASE models. The solid line rep- 
resents the evolution of a galaxy which formed all its stars in an initial, instantaneous burst. The dashed line corresponds 
to a galaxy with an exponentially decaying star formation rate, with an e-folding timescale r = 2 Gyr. The vertical line 
represents the division adopted by Barger et al. (1996), at B-R=2; we find this corresponds to D4000w 1.45. 



of the same name defined in § |3.2| . Galaxies with unusually 
strong Wo{HS) indices are defined as either red or blue 
HJ-strong (rHDS or bHDS); these wiU include both K+A 
galaxies and A+em galaxies. Furthermore, normal galax- 
ies in the red half of this plane generally show little or no 
sign of star formation; both the models and real data of 
galaxies redward of D4000=1.45 have Wo{HS)< SA, and 
any contribution from nebular emission should be small. 
Therefore, we lower our Wo{HS) threshold to SA for the 
rHDS galaxies, as only the truncated star formation mod- 
els pass through this part of the plane. 

Both HDS categories will contain A+cm galaxies, which 
have Wo{OII)> 5 A. We discuss the suggestion by Pog- 
gianti et al. (1999), that these are dusty starbursts, in §pt 
whatever their origin, the A-|-em galaxies are not the re- 
sult of truncated star formation, since such activity is still 
present. To isolate those spectra which most likely result 
from truncated star formation, we define two more galaxy 
classes, subsets of the rHDS and bHDS types: 

• Post-Starburst (PSB) Galaxies: We define 
these to be bHDS galaxies with Wo{OII)< SA. The 
PEGASE models shown in Figure O demonstrate 
that such a spectrum requires the termination of a 
short burst of star formation. Following this event, 
the galaxy spectrum will remain of the PSB type for 
about 300 Myr. Truncated star formation without a 
burst (the solid line) causes the spectrum to evolve 
quickly (^ 100 Myr) through this region, and during 
much of this time emission lines will be present, re- 
ducing Wo{HS) to values less than 5 A. Thus, few (if 
any) bHDS galaxies are expected to result from trun- 
cated star formation without an initial starburst. 



• Post-Star Formation (PSF) Galaxies: These 
are rHDS galaxies with Wo{OII)< 5 A, and their 
spectra can be matched with either a recently termi- 
nated starburst episode, or through the abrupt ces- 
sation of longer term star formation. Galaxies evolv- 
ing along either model track will spend about 300 
Myr as a PSF galaxy; from the D4000 and Wo(Hd) 
indices alone, it is not possible to distinguish between 
these two possibilities. 



3.4. Summary of Definitions 

We have now defined six different types of "unusual" 
spectral types, based on models and the indices of local 
galaxies. All six of these types are rare, and can usually 
be matched by models in which star formation is trun- 
cated (T), sometimes following a short starburst (TSB); 
the galaxies with emission lines are not matched by either 
of these models, and are of uncertain origin (U). These six 
types, their definitions and interpretations are summarized 
in Table g 

4. RESULTS 

4.1. Spectral Index Dependence on Cluster -Centric 

Radius 

The measured values of 1^ (i?(5) , D4000 and WJOII) 
in the maximal sample (|: 2.3.3) are shown in Figure n2, as a 
function of cluster-centric radius, i?/i?200- Field galaxies 
are shown in the right hand panel of each graph, plot- 
ted against an arbitrary abscissa. The solid lines in each 
graph are the weighted median value for the field (right 
panel) and every 150 cluster points (left panel), with 2a 
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Fig. 11. — PEGASE models in the D4000-T4^o (.ff (5) plane, for four different star formation histories which are discussed 
in the Appendix. Nebular emission contribution, not included in the models, will lower the Wo{H6) index of any track 
with ongoing star formation, by at least lA. A galaxy with a constant star formation is at the locus indicated by the 
star, an exponentially declining SFR gives rise to the evolutionary track traced by the long-dashed line. Truncating star 
formation in a galaxy with constant star formation at 4 Gyr causes it to evolve along the solid line. The dashed line shows 
an initial, instantaneous burst of star formation, while the dotted line represents a 200 Myr burst, involving 30% of the 
stellar population of a constantly star forming galaxy (all star formation is truncated following the burst). The plane is 
divided into SSB (short starburst), rHDS (red H(5-strong), blue HDS (blue H(5-strong), SF (star-forming), and Passive. 
The local data arc from Kennicutt (1992b, open symbols) and Kinney et al. (1996, filled symbols). 
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Table 2 




Type 


Wo{OII) 

(A) 


Wc{HS) 

(A) 


D4000 


Interpretation 


K+A 
A+em 
bHDS 
rHDS 
PSB 
PSF 


< 5 
> 5 

<5 

< 5 


>5 
>5 
>5 
>3 
>5 
>3 


< 1.45 

> 1.45 

< 1.45 

> 1.45 


T,TSB 

U 

TSB,U 

T,TSB,U 

TSB 

TSB,T 
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Note. — Six spectral type classifications describing rare 
galaxy types which may be the result of recently terminated 
star formation. Column 5 lists compatible interpretations 
of each spectral type: T: Truncated star formation without 
an initial starburst; TSB: Truncated star formation following 
a burst; and U: Uncertain (due to the presence of emission 
lines) . 



jackknife errors (Efron 1981, Efron fc Tibshirani 1986| ). 
In the central regions, all three indices are inconsistent 
with the corresponding field value with at least 2a con- 
fidence, and they become increasingly field-like with in- 
creasing radius. The Wo{HS) index shows a small increase 
from the center of the cluster out to i?200j significant at 
about the 2a level; the decrease in D4000 over this range 
is measurable with about the same confidence. The me- 
dian value of 14^0(0//) is roughly constant throughout the 
cluster, within 2a limits, but clearly the total distribution 
changes more strongly; particularly striking is the near 
total absence of galaxies with Wo{OII)> 20 A within 0.1 
i?200, as shown previously in Balogh et al. (1997, 1998). 
These trends are consistent w ith those observed in th e in- 
dividual clusters Abell 2390 (Abraham ct al. 199q ) and 
MS1621. 5-1-2640 ( [Morris et al. 19981 ), and are qualita- 
tively consistent with a model in which most cluster galax- 
ies evolve passively, and the mean age of the galaxy stel- 
lar population decreases with increasing distance from the 
cluster center. 

4.2. Galaxies in the Wo{OIiyWo{H5)Plane 

We plot Wo{OII) and Wo{H5) for the cluster and field 
(maximal) samples separately, in Figure O. Each panel 
is divided int o t he same regions shown in Figure S and 
discussed in §3.2. The weighted percentage of galaxies in 
each region is shown, and the error associated with each 
number is determined assuming Poisson statistics. This 
error does not account for the uncertainty due to u nequa l 



scatter in this figure; we will address this issue in §4.2.1 



The fraction of galaxies undergoing some type of star- 
formation (whether SF, SSB or A+em type) is 57 ± 4% 
in the field and only 29 ± 2% in the clusters. Despite the 
established B-0 effect, the mean star formation rate in 
these clusters is less than that in the field, as discussed in 
Balogh et al. (1997; 1998). There is an apparent excess 
of K+A galaxies in the clusters compared with the field, 
significant at almost the 3(T level. As we will show in the 
next subsection, most or all of this difference can be at- 
tributed to the scatter in this figure, and to the slightly 
different absolute magnitude distributions of the cluster 



and field samples. In particular, the systematic effects of 
scatter result in an overestimation of the K+A fraction; 
thus, we interpret the cluster fraction of 4.4 ± 0.7% as an 
upper limit. This is consistent with the results of other 
similarly defined samples at these redshifts, as we show in 



§1 



4.2.1. Determining the True K+A fraction 



The large errors on both indices plotted in Figure O 
spread the data out in this plane; this results in non- 
uniform scatter since galaxies from the most populated 
regions (Passive and SF) will preferentially scatter into 
neighboring areas, biasing the fractions in these regions to 
larger values. In particular, the fraction of K+A galaxies, 
which occupy the high-end tail of the Wo{H5) distribu- 
tion, will be overestimated, and we can expect this effect 
to be larger in the cluster sample, where the Passive galaxy 
population is more dominant than in the field. We can 
demonstrate this by assuming that the true distribution 
of Passive galaxies is the same in both the cluster and the 
field though the absolute number of such galaxies in each 
environment is different. Thus, the Wo{H5) distribution 
of galaxies with Wa{OlT)< 5 A should be the same if the 
fraction of K+A galaxies is also identical. In Figure O, we 
show these distributions normalized to unit area in the top 
(field) and middle (cluster) panels. Subtracting the cluster 
bins from the field bins yields the distribution shown in the 
lower panel. This allows easy comparison of the shapes of 
the Wo{H5) distribution in the cluster and field samples: 
a relative excess of cluster galaxies with a given Wo {H6) 
would result in a negative value in the bottom panel of 
this figure. Instead, we see that the fraction of galaxies 
with Wo{H5)> 5A, corresponding to the K+A galaxies, 
is not significantly different in the two samples. That is, 
the difference between the (normalized) cluster and field 
samples in the two bins with Wo{H5)>h A is equal to zero 
within the la errors. There is a la excess of field galaxies 
with Wo{H5)Ki 3 A, which suggests that even "passive" 
galaxies in the field may have somewhat younger stellar 
populations than similar galaxies in clusters. 
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Fig. 12. — Spectral indices plotted as a function of cluster-centric radius. The sample displayed in the right hand panels 
is the field galaxy sample, plotted against an arbitrary abscissa for display purposes. Positive values of Wo{OII) and 
negative values of Wo{H6) represent emission in these indices. The solid lines are the median index value for every 150 
cluster galaxies, and the 2a jackknife error estimates. The range of data shown is restricted to clarify the characteristics 
of the bulk of the data; arrows indicate data which lie outside the plotted regions. 
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Fig. 14. — The distributions of Wo{HS) for the subsample of galaxies with Wo{OII)< 5A are shown for the field {top 
panel) and cluster {middle panel), normalized to unit area. The number of galaxies in each sample is shown in the top 
right corner; all error bars are la, assuming Poisson statistics. The difference between these two distributions is shown 
in the bottom panel. There is no excess of cluster galaxies with Wo{HS)> 5A, which correspond to the K+A galaxies, 
relative to the field. 



We will now estimate a correction for the fractions in 
each of the regions of Figure |l^. The only way that a pre- 
cise correction can be made is if the a prior distribution 
is known, which is not the case here; therefore, some as- 
sumptions about this distribution will have to be made. 
We will consider two different methods of obtaining this 
correction: the first is applicable to each galaxy class, but 
requires stronger assumptions, while the second method 
provides a more robust estimate but is only valid for the 
K-f-A galaxies. 

Method 1: Monte Carlo Simulations: We could 
estimate the effects of scatter by adding appropriate noise 
to a "pure" sample with no scatter. We attempt to ap- 
proximate such a sample by selecting galaxies with rea- 
sonably high signal-to-noise measurements (S/N> 5) and 
errors in both indices which are below the 25th percentile 
{AWo{OII)< 2.5A and AWo{H6)< 1.5k). Admittedly, 
this is not quite what we want, as it is quite possible that 
this subsample (generally with higher luminosities) is not 
distributed like the full sample; however, it should give 
an approximately accurate correction, as confirmed by the 
second method described below. These data are shown in 



the bottom panel of Figures |15| and |16| for the cluster and 
field samples, respectively. 

We now simulate noisy cluster and field samples by 
adding Gaussian random noise, with a variance chosen ran- 
domly from the la uncertainties of the full (cluster or field, 
as appropriate) data sample. We then generate several re- 
alizations, each one containing a number of points equal 
to the number in the true cluster or field sample. One 
such simulation is shown in the top panels of Figures Hq 
and |l6|. From twenty realizations, we determine, for each 



galaxy class, the average difference between the fraction of 
such galaxies in the noisy (simulated) and true (low-error 
sample) distributions, and the variance of that difference. 
We will refer to this value as the "scatter correction" . 

In Table ^, we show the measured fraction of each galaxy 
type, from Figure O, in columns 2 (cluster) and 5 (field). 
Columns 3 and 6 show the estimated scatter correction. 
Subtracting these numbers from the appropriate measured 
value, including the error in this correction, yields the cor- 
rected fractions in columns 4 and 7. 

We note that, in Table |[ the corrected K-l-A fraction is 
equal to zero within 2a uncertainties. However, this does 
not mean that K-f A galaxies are not present in the sam- 
ple; to provide concrete examples we show, in Figure [17| , 
sample spectra for two K-l-A galaxies: one found in the 
field and one in the cluster. Note that the Ca II K line, 
prominent in the bottom spectrum, is weak/absent in the 
top spectrum; this is an example of two galaxies which 
are assigned the same classification (K-f A), though their 
spectra are not identical. 

Method 2: Gaussian Cloud Analysis: An al- 
ternative way to determine the overestimate of the K-|-A 
fraction is to assume that it is due only to scatter in the 
Wo{HS) direction from the Passive galaxies, i.e. those 
with Wo{H6)< SA and Wo{OII)< 5 A. We can represent 
the Wo {H5) index of each galaxy as a Gaussian probabil- 
ity function, with variance given by the la index uncer- 
tainty. For each galaxy in the Passive sample, we calculate 
the probability that Wo{HS)> 5 A, and then compute the 
weighted sum of all these probabilities. A smaller amount 
of "backscatter" (from the K-l-A region to the Passive) 
will also take place, and this must be subtracted off in an 
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Fig. 15. — The bottom panel shows cluster data with S/N> 5, and errors in Wo{H6) and Wo{OII) less than 1.5A and 
2.5A, respectively. Error bars are la. We add noise, characteristic of the full sample error distribution, to each point, and 
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shown in the top panel. The plane is divided into the same regions defined in Figure 0; the numbers in each panel are the 
weighted percentage of objects in that bin. This experiment is used to calculate the scatter correction given in Table Q. 
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Fig. 16. — This is the same as Figure |l^, but for the field sample. Again, the results are tabulated in Table ||. 
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Table 3 
Systematic Error Estimations in the Wo{OII)-Wo{H5) Plane 



Object Type 




Cluster (%) 






Field (%) 






Measured 


Scatter 


Corrected 


Measured 


Scatter 


Corrected 


K+A 


4.4 ±0.7 


2.8 ±0.5 


1.6 ±0.9 


1.9 ±0.6 


2.2 ±0.9 


-0.3 ±1.1 


A+cm 


4.4 ±0.7 


1.4 ±0.5 


3.0 ±0.9 


10.3 ±1.7 


4.0 ±1.3 


6.3 ±2.1 


SF 


18.1 ±1.6 


-2.4 ±0.6 


20.5 ±1.7 


39.0 ±3.3 


-3.3±2.2 


42.3 ±4.0 


Passive 


67.0 ±2.6 


-7.1 ±1.0 


74.1 ±2.8 


40.7 ±3.0 


-6.2±1.4 


46.9 ±3.3 


SSB 


6.1 ±0.8 


5.0 ±0.8 


1.1 ± 1.1 


8.1 ±1.5 


3.4 ±1.4 


4.7 ±2.1 
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Fig. 17. — Two high S/N examples of galaxies in our sample with a K+A type spectrum, in the rest frame. The bottom 
panel shows a cluster member in Abell 2390 at z = 0.2211 (S/N=22, ppp ^^100537); the top panel shows a field galaxy 
in the field of MS1512 at z = 0.4926 (S/N=7, ppp #200607). These are good quality spectra, and not characteristic of 
typical spectra in the CNOCl sample. Note that, although both galaxies are classified as K+A, the top spectrum lacks 
the strong Ca II K line present in the bottom spectrum. 
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analogous manner. This difference would give the total 
fraction of galaxies which have scattered from the Passive 
region into the K+A region, if the distribution of data 
in this plane was representative of the error-free distri- 
bution; since the uncertainties significantly blur this dis- 
tribution, this difference results instead in the maximum 
amount of K-f A "contamination" expected. Following this 
procedure, we determine that the maximum correction to 
the K-l-A fraction should be 2.8% for the cluster, and 2.5% 
f or the field (appropriate for the maximal data sample of 

We can now estimate the minim,um correction required 
by counting the weighted number of Passive galaxies with 
Wo{H5)< — SA. Since the mean Wo{H5) in the Passive 
region is greater than zero, and we expect few if any Pas- 
sive galaxies to actually have Wo{H6)< —5, this fraction 
of galaxies represents the minimum amount of scatter into 
the K-l-A region (without a backscattcr correction). For 
the maximal cluster sample, this minimum correction is 
1.9%; for the field, it is 1.2%. 

If we adopt the average of the minimum and maximum 
corrections, and half the difference to represent (arbitrar- 
ily) a 2cr uncertainty, we determine the best correction 
to be 2.3±0.25% for the cluster, and 1.8±0.3% for the 
field, resulting in corrected K-l-A fractions of 2.1±0.7% 
and 0.1±0.7%, respectively. These corrections are some- 
what smaller than the ones in Table 0, but equivalent 
within the uncertainties. Although we take this method 
to be the more robust estimate for the K-l-A correction, it 
confirms that the corrections listed in Table || are likely to 
be reasonable estimates. 

Luminosity Limit: In Figure |8|we show the M^ dis- 
tribution of the K-|-A sample, compared with the overall 
distribution of the maximal sample. This figure shows that 
many of the K-l-A galaxies lie below Mr = —18.8-1-5 log h, 
where the magnitude weights are large, and the maximal 
sample is not statistically representative of a complete pho- 
tometric sample. Thus, any absolute evaluation of the 
K-l-A fraction will be sensitive to the magnitude limit and 
statistical correction. In consideration of this, wc show the 



results for the luminosity limited sample (§2.3.2) in Figure 
19, For this restricted sample, the fraction of K-f A galax- 
ies is equal to 3.3 ± 0.7% in the cluster, and 2.7 ± 0.8% in 
the field; the two are equivalent within Icr uncertainties. 
Adopting the Gaussian cloud scatter correction (Method 
2), we find that these two numbers are overestimated by 
1.8±0.4% and 1.5±0.3%, respectively; thus, the true frac- 
tion is reduced to 1.5±0.8% in the cluster, and 1.2±0.8% 
in the field. 

Summary of Corrected Values: Although the 
amount of scatter correction necessary cannot be pre- 
cisely determined (since the true distribution of the data 
is unknown), the uncorrected, cluster K-f A fraction of 
4.4±0.7% is a secure upper limit. Our best estimate of 
the true K-fA fraction is thus obtained by applying the 
corrections determined by the Gaussian cloud estimate. 
For the maximal galaxy sample, we find the fraction of 
K-fA galaxies is 2.1±0.7% in the cluster, and 0.1±0.7% in 
the field. For the luminosity limited sample, the fraction 
is 1.5±0.8% (cluster) and 1.2±0.8% (field). We conclude 
that, once scatter and luminosity effects are accounted for, 
there is no significant difference between the frequency of 



K-f A galaxies in the cluster and field environments. We 
will compare these numbers with the results of other sur- 
veys, both locally and at similar redshifts, in §|[ 

4.2.2. The Radial Dependence 

The fraction of each galaxy type is plotted as a function 
of distance from the center of the cluster, normalized to 
i?200 5 in Figure |20|. The three "active" classes SF, A-fem 
and SSB have been grouped together for clarity. We have 
corrected each value for the estimated effects of scatter, us- 



ing Method 1 in §4.2.1, and included the uncertainty of this 



correction in the error bars. The galaxy population within 
the cluster exhibits a steadily increasing fraction of pas- 
sively evolving galaxies and a steadily decreasing fraction 
of star forming galaxies with decreasing radius. This is at 
least partly expected from the morphology-density rela- 
tion (Dressier 1980) but, as shown in Balogh et al. (1998) 
this relation cannot completely account for the decrease 
in star formation rate toward the center of the cluster. 
The fraction of K-f A galaxies within the cluster is never 
significantly in excess of the field value, at any radius. 

4.2.3. The Redshift Dependence 

The fifteen clusters in the CNOCl sample cover a large 
range in redshift, from 0.18 to 0.55 and, thus, it is of inter- 
est to investigate the possible redshift dependence of the 
populations. To ensure an equitable c ompar ison, we con- 
sider the luminosity limited sample (§ |2.3.2 ) here (recall 
that an evolution correction is applied to the luminosi- 
ties). First, we show in Figure |2| the redshift distribution 
of this sample, separately for the cluster and field. This 
figure shows that the cluster observations are well divided 
into three, well populated redshift bins: 0.15 < z < 0.28; 
0.28 < z < 0.35; and 0.35 < z < 0.7. 

In Figure E3we plot the fraction of each galaxy type, 
as defined in § |3.2| and Figure |l^, in each of the three red- 
shift bins described above. Again, we have grouped all 
star-forming galaxies together (SF, A-fem and SSB) for 
clarity, and corrected the val ues for the systematic effects 
of scatter using Method 1 of §4.2.1. The cluster sample is 
shown in the bottom panel, and the field in the top panel. 
A small increase in the fraction of star forming galaxies 
with redshift is observed, though only significant at about 
the Icr level. Importantly, the abundance of K-f A galaxies 
does not show an increase with redshift, which shows that 
the smallness of the fraction we observe is not just due 
to the fact that our sample is dominated by low redshift 
z w 0.2 galaxies. 



4.2.4. The Near-Field Galaxies 

In Figure [2^ we show the Wo{OII)-Wo{H5) relation for 
t he ne ar-field galaxies in the luminosity limited sample 
(§ |2.3.2| ). The near-field galaxies are a subset of galaxies 
with velocities intermediate between our cluster and field 
definitions (between 3 and 6 times the cluster velocity dis- 
persion), and will include both galaxies that are far from 
the cluster (« 50 Mpc), with velocity differences that re- 
flect the Hubble flow, and infalling galaxies with peculiar 
velocities in excess of three times the cluster velocity dis- 
persion. None of the galaxy abundances are significantly 
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Fig. 18. — The k- and evolution-corrected absolute magnitude distribution for the K+A galaxies {solid line) and the 
full maximal sample {dashed line). Many of the K+A galaxies are less luminous than Mr = —18.8 + 5 log ft., where the 
statistical corrections are large and the maximal sample is not representative of a complete photometric sample; the 
abundance of K+A galaxies is therefore sensitive to the magnitude limit and statistical correction used. 
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Fig. 19. — The same as Figure O, but for the luminosity limited sample described in §2.3.2. The fraction of K+A galaxies 
in the cluster and field samples are equivalent within la uncertainties. 
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Fig. 20. — The weig hted fraction of each type of galaxy, as selected based on its position in the Wo{OII)-Wo{H5) plane 
shown in Figure O, as a function of project ed rad ius. The fractions and la error bars are corrected for the systematic 



effects of nonuniform scatter, as discussed in §4.2.1. The connected points represent the cluster sample; the isolated points 



at i?/i?200 — 3.6 represent the values in the field sample, plotted at an arbitrary radius for display purposes only. 
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Fig. 21. — The redshift distribution of the cluster (solid line) and field (dashed line) samples, limited to galaxies more 
luminous than Mr ~ —18.8 + 5 log ft, (corrected for band shifting and evolution). 
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Fig. 22. — This figure shows the fraction of each galaxy population defined in § |3.2| and Figure |1^, as a function of redshift. 
Only galaxies more luminous than Air = —18.8 + 5 log ft, are considered. The field galaxy sample is shown in the top 
panel, and the cluster sample in the bottom panel. The stati stical ly weighted data are corrected for the systematic effects 
of scatter in the Wo{H6)-Wo{OII) plane, as discussed in § 4.2.1 . All error bars are la. The K+A fraction is low in all 
three redshift bins. 



different from the field (c.f. Figure |I^. Thus, even consid- 
ering galaxies which may be in the infall regions of clusters, 
we do not detect evidence for an excess of star forming or 
K+A galaxies, relative to the field. 



dant, relative to the FSB galaxies, by a factor of ~ 5 — 10. 
Though this number is quite uncertain, if the large excess 
is real it may have significant consequences, as we discuss 
in §1 



4.3. HDS and PSF Fractions 

The measu remen ts of Wo (HS) and D4000 for the maxi- 
mal sample (§ 2.3.3| ) are shown in Figure p4, for the cluster 
and field, separately. Each plane is divided into the five re- 
gions shown in Figure ni| the number in each panel is the 
weighted percentage of galaxies within it. The numbers 
in brackets, for the cluster sample, are from Barger et al. 
(1996), corrected for small difference s in region definitions. 
We will compare these numbers in §5.2. 



We can correct the fraction of galaxies in each region of 
this figure for the effects of scatter in a manne r analogous 
to the Monte Carlo method described in §4.2.1; the results 
are shown in Table 0. After this correction, both bHDS 
(at 1.0±1.2%) and rHDS galaxies (3.6±1.4%) in clusters 
are not more common than they are in the field (2.7±1.8% 
and 9.4±2.2%, respectively). 

As discussed in §3.3, the HDS galaxies include emission 
line galaxies; populations with recently truncated star for- 
mation should fall into either the FSB or PSF subclasses. 
In Table | we show, for the cluster and field separately, the 
"raw" PSB and PSF measurement, the estimated scatter 
correction (assuming that the relative scattered fraction is 
the same for the PSF/PSB classes as for the rHDS/bHDS 
classes) and the corrected estimate of the "true" fractions. 
Again, the fractions of PSB and PSF galaxies are not sig- 
nificantly more common in clusters than they are in the 
field. In both cases, the PSF galaxies are more abun- 



4.3.1. Reddening, Metallicity and IMF Effects 

A striking feature of Figure g4| is the abundance of galax- 
ies, both cluster and field, with D4000 ^ 1.9 and Wo{Hd) 

;> 2A, which are not matched by any of the models in Fig- 
ure nil This problem has also been pointed out by Couch 
& Sharpies (1987), Morris et al. (1998) and Poggianti & 
Barbaro (1996), among others, and persists if photometric 
colors are considered instead of D4000. The problem has 
been substantially alleviated by considering the narrower 
D4000 index but it still persists. 

In Figure ES, we plot only the data with S/N >15 and 
uncontaminated by bright night sky lines near the Wo {US) 
and Wo{OII) indices. For reference, we plot the initial 
burst, Salpeter IMF model shown in Figure n^, with a 
solid line. Using the deredden task within IRAF, we red- 
den this model by a large amount, Ej^_y—0.5; this moves 
the model to the right by about 0.15, as represented by the 
dotted line. This reddened model extends in D4000 nearly 
to the full extent of the data, and the scatter in Wo {HS) 
at the red end is fairly uniform about it. 

We also consider two other model variations. The short- 
dashed line in Figure pa represents an initial burst model 
for a galaxy with super-solar metallicity, Z=2.5.^o, gener- 
ated from the updated Bruzual & Chariot (1993) mod- 
els. Although this model does indeed extend to very 
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Fig. 23. — Selected from the luminosity limited sample, we show the near-field galaxy population, which has velocities 
intermediate between our cluster and field definitions, in the Wo{OII)-Wo{HS) plane, with la error bars. The regions 
are the same as in Figure |9|, and the number shown in each region represents the weighted percentage of galaxies within 
that region. The error is estimated assuming Poisson statistics. The fraction of galaxies in each class is not significantly 
different from the field values shown in Figure p^. 



Table 4 
Systematic Error Estimations in the D4000-Wo{HS) Plane 



Object Type 




Cluster (%) 






Field (%) 






Measured 


Scatter 


Corrected 


Measured 


Scatter 


Corrected 


bHDS 


4.7 ±0.8 


3.7 ±0.9 


1.0 ±1.2 


7.2 ±1.3 


4.5 ±1.3 


2.7±1.8 


rHDS 


9.7±1.0 


6.1 ±1.0 


3.6 ±1.4 


11.1 ± 1.8 


1.7 ±1.3 


9.4 ±2.2 


SF 


10.5 ±1.1 


-12.1 ±1.1 


22.6 ±1.5 


32.7 ±3.0 


-11.3 ±2.2 


44.0 ±3.7 


Passive 


73.5 ±2.8 


-0.8±1.2 


74.3 ±3.0 


44.7 ±3.1 


1.7 ±1.6 


43.0 ±3.5 


SSB 


1.5 ±0.4 


3.1±0.7 


-1.6 ±0.8 


6.3 ±1.5 


3.2 ±1.4 


3.1 ±2.0 



Table 5 
Abundances of PSB and PSF Galaxies 



Object Type 


Measured 


Cluster (%) 
Scatter 


Corrected 


Measured 


Field {%) 
Scatter 


Corrected 


PSB 
PSF 


1.3 ±0.4 
4.9 ±0.6 


1.0±0.2 
3.0±0.5 


0.3 ±0.4 
1.9 ±0.8 


0.7 ±0.3 
3.6±0.9 


0.4±0.1 
0.5±0.5 


0.3 ±0.3 
3.1 ±1.0 
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Fig. 24. — The data presented in the T) 4000- Wo (HS) plane, divided into a field and cluster sample. The models and region 
delineations over-plotted are the same ones shown in Figure O. The number in each region is the weighted percentage 
of galaxies in that region, and the errors assume Poisson statistics. The numbers in brackets, for the cluster sample, are 
from Barger et al. (1996), corrected for the small differences in region definitions. Our results are compared with those 
of Barger et al. in §5.2. The sample error bars represent the mean la error. 
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large D4000, the Wo{H6) index is always negative^^ for 
D4000> 1.9 and, thus, lies below the bulk of the data. Sec- 
ondly, we consider the effects of an extreme IMF model, af- 
ter Rieke et al. (1993), in which no stars less massive than 
2Mq are formed. An initial burst of star formation with 
this IMF is shown in Figure ^ as the long-dashed line. 
This model is also able to match the strongest Wo{H6) 
indices of our data with D4000> 1.9. However, the stellar 
population resulting from such a burst is very short lived, 
about 300 Myr; thus, it is unlikely that many of these red 
galaxies can be undergoing such a burst, and we expect 
that significant dust reddening is the best explanation of 
the data. This supports the finding of Poggianti et al. 
(1999), that dust obscuration in galaxies at these redshifts 
plays an important role in the appearance of their spectra. 



5. COMPARISON WITH PREVIOUS WORK 
5.1. The Low Redshift Universe 

The Las Campanas Redshift Survey (LCRS), as pre- 
sented in Zabludoff et al. (1996), is a large (- 24,000 
galaxies), unbiased sample of relatively nearby galaxies 
(0.05 < z < 0.13) for which good quality spectra and line 
indices are available. Zabludoff et al. identify 21 of their 
11,113 galaxy subsample (Mr > -19 -I- blogh) as E-l-A 
(renamed here K-t-A) types, based on Wo{OII) and an 
average equivalent width of US, II7, and II/3 (< HSj(3 >). 
They estimate that 1785 galaxies in their sample lie within 
the infall radius of a rich cluster, and five of these are iden- 
tified as K-l-A types; therefore, 0.28 ±0.12% of the cluster 
population are K-l-A galaxies, compared with 0.17±0.04% 
in the field. These fractions are equivalent within la un- 
certainties. The definition of a K-l-A galaxy adopted by 
Zabludoff et al. differs slightly from the one we have 
adopted here (§^^) in two respects. First, they use a min- 
imum Balmer index of < H5^(3 >=5.5A. From the PE- 
GASE models we find that, for Wo{H5)k 5A, < HS-^P > 
-Wo(HS)^ 0.7; thus, our limit of Wo{H6)>5 A is nearly 
equivalent to theirs. Secondly, Zabludoff et al. require 
Wo{OII)< 2.5A; our lower S/N requires us to use a more 
generous definition of Wo{OII)< 5A. Adjusting for this 
difference allows the inclusion of an additional ^ 13 K+A 
galaxies in the LCRS sample, raising the total K-l-A frac- 
tion to 0.30 ±0.05%. Since the uncertainties on the LCRS 
line indices are small (due to the high S/N o f the data), 
a scatter correction like the one discussed in § [4.2.1 is not 
made. 

In our luminosity limited sample (§ 2.3.2| ), which most 
closely matches the luminosity range of the Zabludoff et al. 
sample, the "raw" K±A fraction in the field is 2.7±0.8%, 
which we regard as an upper limit. Our best estimate 
of this fraction, corrected for scatter using the Gaussian 
cloud method (§ [4.2.1| ), is 1.2±0.8%. Thus, the evolution 
in the K±A fraction over this fairly small redshift range 
(from z ~ 0.1 to z ~ 0.3) amounts to a factor of only 
about 4±2.7, and the two fractions are actually consistent 
at about the la level. 



5.2. Couch & Sharpies (1987) and Related Work 



Much of the original work in this field was done by 
Couch & Sharpies (1987), whose dataset of 152 galaxies 
more luminous than about Mr = —20 + 51og/i in three 
clusters at z « 0.31 was later analyzed in more detail by 
Barger et al. (1996). As shown in Figure g^, we find some- 
what fewer rHDS and bHDS galaxies than Barger et al., 
though the difference is significant at less than the 2a level, 
and may be partly due to an impreci se m apping of B-R 
to D4000 (which was determined in |3.3| ). Barger et al. 
performed a useful and detailed analysis of the distribu- 
tion of galaxies in this figure, and compared it with model 
simulations to determine the importance of starbursts to 
cluster galaxy evolution. They suggest the data are well 
represented by a model in which ~ 30% of cluster galax- 
ies have undergone a 0.1 Gyr starburst (which produces 
10-20% of the the galaxy's final stellar mass) in the last 
~2 Gyr. Given the small number of galaxies and the con- 
siderable uncertainties on the line indices, these numbers 
are uncertain, as the authors acknowledge. In §0 we will 
argue from the fraction of K+A galaxies in the CNOCl 
sample that the importance of such short-lived starbursts 
in clusters is unlikely to be this large. 

The Couch & Sharpies data are also considered by Pog- 
gianti & Barbaro (1996), who conclude from their model- 
ing that strong starbursts are required to match the K±A 
galaxies in that sample. However, the comparison between 
the data and models may be compromised by a differ- 
ence in index definition. The Couch & Sharpies data have 
Wo{H5) indices as large as 9 A, with many points around 
5 <Wo{H5)< 6 A, while the post-starburst models of Pog- 
gianti & Barbaro struggle to reach Wo{H5)= 6 A (see 
their fig. 4), and then only for short periods of time. By 
dividing the data and models into broad categories (i.e. 
Wo{H5)> 3 and {B — R) < 2) it is possible to misinter- 
pret the results if the models do not match the data within 
a given category, as may be the case here. 

5.3. C11358+62 and CFRS 

Fisher et al. (1998) analyzed spectra for 232 galaxies 
within a ~4 Mpc region around the z « 0.32 cluster CI 
1358+6245, with a magnitude Umit of Mr w -19.1 + 
5 log h. This improves on the preliminary work presented 
by Fabricant et al. (1991), which used only 70 spectra 
in the central regions of this cluster. These authors de- 
fine K+A (which they term E+A) galaxies in a manner 
similar to Zabludoff et al. (1996), but with more gen- 
erous limits in < HSjf3 > and 11^0(0//). They find 11 
(4.7±1.9%) of their cluster galaxies show K+A type spec- 
tra, remarkably similar to our total cluster fraction, before 
the necessary scatter correction. Since their uncertainties 
are also comparable to ours, we expect this correction also 
to be similar, resulting in a reduced, true K+A fraction. 

A similar comparison can be made with the work of 
Hammer et al. (1997), who find a K+A fraction of 
4.9±1.5% in the CFRS field galaxy sample, with a less 
precise K+A definition, roughly similar to that of Fisher 
et al. Again, this number should be reduced due to the 
overestimation resulting from unequal scatter and, thus, 
the amount of evolution implied relative to the low red- 
shift (LCRS) sample will be significantly less than the or- 
der of magnitude implied by the raw fraction alone. This 



^^Recall, from §2.2, that a negative Wo{H5) index does not necessarily imply an emission feature. 
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Fig. 25. — PEGASE star formation models presented in the D4000~W^o (-ff <5) plane, with regions delineated as in Figure 
11 The solid line is the standard, Salpeter IMF initial burst model, also shown in Figure |ll|. This model, reddened by 
i?5_y=0.5, is shown as the dotted line. A truncated IMF initial burst (Rieke et al. 1993) is shown as the long-dashed 
line, and a high metallicity (GISSEL96) model is the short-dashed line. The data plotted are only those with high 
signal-to-noise ratios (S/N>15) and uncontaminated by bright night sky lines; the error bars are la. 
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number is similar to our raw cluster K+A fraction, but 
somewhat larger than our field value. Thus, it supports 
our conclusion that K+A galaxies are not preferentially 
found in clusters; however, we note that the mean redshift 
of the CFRS, z k. 0.6, is considerably higher than that 
of the CNOCl sample, and evolutionary effects may be 
important. 

5.4. The MORPHS Collaboration 

After this paper was submitted, the MORPHS collab- 
oration published a similar study based on ten clusters 
at 0.37 < z < 0.56, including 657 galaxies brighter than 
Mr ~ -19 + 51og/i. (Dressier et al. 1999; Poggianti et 
al. 1999). These authors conclude that galaxies without 
detectable [O II] emission but Wo{H5)> 3 A (which they 
call k+a/a+k) are significantly more common in the clus- 
ter sample than the field, (21±2% compared with 6±3%). 
Both the largeness of these fractions, and the significance 
of their difference, are apparently in conflict with the re- 
sults we have presented in ^.2| . This difference warrants 
a detailed comparison of the data and analysis, which we 
present in this subsection. 

We will first show that our data are of comparable qual- 
ity to the MORPHS data, and then suggest that the dif- 
ferences in our conclusions may result from (1) different 
selection criteria (both of galaxies and of clusters); (2) dif- 
ferent treatment of uncertainties; and (3) differences in 
galaxy classifications and definitions. 

5.4.1. A Comparison of Data Quality 

We have obtained the publ icly available MORPHS data 
from http://www.ociw/'-^irs, and computed our indices 
from their data directly. However, the MORPHS data 
do not include error vectors, and uncertainty estimates of 
their indices are unpublished; thus care must be taken be- 
fore quantitative comparisons are made with our results. 
In particular, continuum S/N ratio s ca nnot be computed 
in the same way as we have done in j |2.1| . Instead, we define 
an alternate estimate as the mean flux per pixel divided by 
the r.m.s. in the wavelength region 4050 < A/A< 4250. 
This underestimates the true S/N, especially for high S/N 
spectra, due to real features in the spectrum. In Figure 
we show how this "alternate" value of S/N, computed 



made on the MORPHS sample by correlating the errors on 
the CNOCl data with our alternative S/N measurement. 
From the S/N measured on the MORPHS data, we can 
then determine approximate index uncertainties. Repre- 
senting the uncertainties in Wo (HS) and Wo (Oil) by AHS 
and A[0 II], respectively, we find AHS = 9.2 x {S/N)-^-°^ 
and A[0 II]= 32.8 x {S/N)-^-^^. 

The uncertainty of a line index will be reduced if the 
resolution is improved, and the continuum S/N ratio is 
fixed. The MORPHS spectra were obtained using differ- 
ent telescopes and instruments, and are not all of the same 
resolution; the best resolution spectra (obtained with the 
COSMIC spectrograph) have a dispersion of 3.1 A per 
pixel, and are smoothed to the instrumental resolution of 
~ 8 A, and rebinned to 10 A per pixel. Thus, the aver- 
age FWHM of emission lines is comparable to that of the 
CNOCl spectra, ~ 16.5 A, as confirmed by the measure- 
ment of several [O II] lines. 

A second property which must be tested is the equiv- 
alence of the CNOCl line indices to those measured by 
the MORPHS collaboration. In Figure 128 we show the 
values of our Wo {HS) index, defined in §[2.2| and measured 
on the MORPHS spectra, compared with those of Dressier 
et al. (1999). We exclude galaxies with quality indices of 
4, or flagged as uncertain. The latter were usually com- 
puted by interactively fitting the continuum and then fit- 
ting a Gaussian profile to the absorption or emission lines. 
Given the noisy nature of the data, and the generally com- 
plex region of the spectrum under consideration, we find 
that such an interactive approach is difficult to duplicate, 
and makes uncertainty estimates difhcult to determine re- 
liably. Nonetheless, the correlation between our HS index 
and that of Dressier et al. is fairly good. In particular, we 
recover the large Wo (HS) indices measured in many cases 
by Dressier et al., showing that our index is not insensitive 
to spectra of this type. The scatter in this figure reflects 
the fact that both of these indices have considerable un- 
certainties. 



5.4.2. Reanalysis of MORPHS Data 
In Figure |29| we show the MORPHS data in the 



the CNOCl galaxies, compares with our definition in §2.1 
The correlation is evident, though the scatter is large, and 
the underestimation of S/N at high S/N values is clear. 
Based on this correlation, we will assume that our alter- 
native definition can hold as a reasonable proxy for the 
more realistic measurement used to qualify the CNOCl 
data. 

In Figure E^, we compare the S/N distribution of the 
CNOCl sample to that of the MORPHS sample (exclud- 
ing from the latter galaxies that have a quality index of 4, 
which are defined as those with S/N sufhcient only for a 
redshift determination). In both cases, S/N is computed 
in the "alternative" manner described above. Both sam- 
ples have similar distributions, apart from a small excess 
of low S/N objects in the MORPHS sample^^ We can 
crudely assign an uncertainty to the index measurements 

^^These S/N ratios arc computed per pixel, and the CNOCl spectra are sampled at 3.45 A per pixel, while many of the MORPHS spectra 
are more coarsely sampled (for example, the COSMIC spectra are rebinned to 10 A per pixel). 



Wo { OII)-Wo (HS) plane, with our region delineations from 
§ 3.2| , and the measurements we have made using our index 
dehnitions (§p^). We exclude data of quality index 4, and 
mark as open circles those galaxies which are either flagged 
as having uncertain Wo{HS) measurements, or for which 
no Wo{HS) measurement is listed in the MORPHS cata- 
logue (i.e., they are listed as either INDEF, or 0). Further- 
more, no magnitude limit is imposed. The sample error 
bars represent the mean error in each index. The field sam- 
ple is somewhat biased to lower S/N galaxies and, hence, 
results in larger mean index errors. 

We will not attempt to evaluate the fractions of each 
type of galaxy (K-l-A, SF etc.) that we have defined for 
galaxies in this plane, due to the difficulty in correcting for 
the morphologically based selection (see Small et al. 1997). 
There is clearly, however, a strong presence of K-l-A galax- 
ies in the cluster sample that is not observed in the field. 
For most (20/24) of the galaxies with Wo{HS)> 10 A, 
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5.4.1. The solid line, drawn for reference, is where the points would like if both S/N 

The alternate 

definition underestimates the S/N, especially at high values, where real features in the continuum add to the r.m.s. which 
is used instead of the noise. 



alternate definition described in 

measurements were equal. Most of the data show a strong correlation, though there is significant scatter 
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Fig. 27. — We compare the S/N distribution for galaxies in the MORPHS sample with quality indices less than A(solid 
line) to that in the CNOCl sample (dashed line). In this case, S/N is measured as the ratio of the mean flux per pixel in 
the range 4050 < X/A < 4250 to the r.m.s. in the same range. The bottom panel shows the number of galaxies in each 
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Dressier et al. (1999) either did not measure Wo{HS) or 
the measurement was flagged as uncertain (open circles). 
Often, this is due to poor sky subtraction near the line, 
which should result in large index uncertainties that justify 
their remo val fro m the sample (as we do for our maximal 
sample in §2.3.3). However a large population of cluster 
K+A galaxies with Wo{HS)< lOA persists. This strong 
population is clearly not as abundant in the MORPHS field 
sample, in contrast with the results found in our sample 
(Figure O). Due to the unusual selection criteria, how- 
ever, the absolute fraction of K+A galaxies in the two 
samples cannot be easily compared in a quantitative way. 



5.4.3. Possible. Smi.rr.p.s of mffe.re.nr.e.s Between MORPHS 



and CNOCl 



Given that the CNOCl resolution, S/N distribution 
and line index definitions are comparable to those of the 
MORPHS collaboration, what is the source for the differ- 
ence in results? We discuss three possibilities below. 

Cluster Selection: Given that galaxy clusters are a 
very heterogeneous class of objects, with differing galaxy 
populations, some or all of the difference between the 
MORPHS and CNOCl results may be in cluster selection. 
The CNOCl clusters are among the most X-ray luminous 
at these redshifts, often containing cooling flows, and are 
mostly dynamically relaxed (Lewis et al. 1999). It is quite 
possible that, as a class, these clusters differ from the more 
heterogeneously selected clusters in the MORPHS sample 
(though we confirm K+A fractions of < 5% for the two 
clusters in common with both samples). Furthermore, the 
MORPHS clusters are generally at higher redshifts than 



the CNOCl targets; all but four of the CNOCl clusters lie 
at redshifts less than that of the lowest redshift MORPHS 
cluster, A370. 

Galaxy Selection: The selection of spectroscopic tar- 
gets from the photometric sample differs in the two sur- 
veys. The MORPHS sample was primarily selected based 
on galaxy morphology; Sd/Irregular galaxies are oversam- 
pled relative to earlier types. Since late-type galaxies were 
targeted preferentially, the detection of a large population 
of starburst and/or post-starburst galaxies may not be 
surprising. Dressier et al. (1999) attempt to correct for 
this by comparing the morphological distribution of their 
spectroscopic sample with the photometric sample, but 
this correction is not completely reliable since the st atisti- 
cal correcti on for late type field galaxies is uncertain (Small 
et al. 1997). Alternatively, the primary selection criterion 
in the CNOCl sample is apparent magnitude, with small 
corrections for galaxy color and position. In particular, we 
note that any color selection effect is quite small, and there 
does not seem to be a significant bias to our results as a 
function of galaxy spectral type at th ese apparent mag- 
nitudes. As discussed in section 2.3.1, we used weighting 



functions to correct our sample of line indices to be rep- 
resentative of the full photometric sample in the cluster 
fields. 

We also note that the CNOCl survey includes galaxies 
at much larger distances from the cluster center, as shown 
in Figure Bft However, the lack of an observable gradient 
in the K-fA fraction (Figure E3) makes it unlikely that 
this can have a strong effect. 

Uncertainties: No uncertainties are published for the 
line indices presented in Dressier et al. (1999) or Poggianti 
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Fig. 29. — Our measurements of Wo{OII) and Wo{H6) on the MORPHS cluster and field samples. Data which are 
of quality index 4 have been excluded. Open circles are data for which either no Wo{H6) measurement is presented in 
the MORPHS catalogue, or the measurement is flagged as uncertain. The sample error bars represent the mean, la 
uncertainty, computed assuming the CNOCl correlation between S/N and index uncertainty holds for the MORPHS 
sample. This figure may be compared with the CNOCl data presented in Figure ^, but note that the selection criteria 
are different in that the MORPHS data are biased toward late type galaxies (which will generally have stronger Wo (Oil) 
and Wo{H6) indices). There is a clear excess of K+A galaxies in the MORPHS clusters, relative to their field sample. 
Although the absolute fraction of K+A galaxies in this sample is difficult to determine reliably, this excess suggests that 
the MORPHS clusters may have different galaxy populations than the X-ray luminous clusters considered in the present 
work. 
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Fig. 30. — The distribution of distances from cluster centers, for gala xies in the CNOCl {solid line) and MORPHS {dashed 



line) samples. The CNOCl distribution is weighted by Wgpoc (see § 2.3.1 ), but renormalized so that the total number of 



galaxies represented in this figure corresponds to the number in the spectroscopic sample. On average, galaxies in the 
CNOCl sample come from radii about 3 times larger than galaxies in the MORPHS sample. 



et al. (1999). As we have shown in § 4.2.1, the large un- 



certainties of our indices result in an overestimate of the 
K+A fraction; since we believe the MORPHS index uncer- 
tainties to be comparable to ours, the same holds true for 
their sample, and some correction for this effect must be 
made. Index measurements for galaxies with undetectable 
or hard to measure lines are not listed in Dressier et al. 
and, thus, the scatter in the measurements that do exist 
does not accurately reflect the average uncertainty of the 
data. 

Galaxy Classifications: Based on unpublished re- 
computations of the Barbaro & Poggianti (1996) models, 
using their new Wo{H5) definition, Poggianti et al. (1999) 
have chosen a threshold of Wo{H6)> 3 A for their classi- 
fication of k-|-a/a-f k galaxies. This is considerably lower 
than we have adopted for our K-l-A definition, and clearly 
leads to the identification of a larger population of "un- 
usual" galaxies. We have justified our limit not only by 
comparison with the PEGASE models, but also by com- 
parison with the local data of Kennicutt et al. (1992a) 
and Kinney et al. (1996). From Figure || alone, it is clear 
that galaxies with Wo{HS)> 3A and no detectable [O II] 
emission are not exclusively K-|-As, but include normal SO 
and early spiral galaxies. Since the Dressier et al. Wo{H6) 
threshold is significantly lo wer than previously ado p ted by 
many other auth o rs (e.g.. prcsslcr fc Gunn 198^, Fabri 
cant 4t al 
19961 



1991 



Za bludo 
199^ 



prcsi 

tfcTi 



aX 199(:, Abraham et al 



Morris et al. 199^ ), we chose not to adopt their 
lower limit. Adopting 3 A instead of SA as our lower limit 



for the classification of K-l-A galaxies would increase the 
fraction of galaxies that satisfy this definition^^ . To some 



extent, of course, the choice of limiting Wo{HS) is quite 
arbitrary and of little consequence as long as it is inter- 
preted as such, and accounted for when comparisons with 
other work is made. 

5.5. Summary 



et al 



With the exception of the MORPHS results (Dressier 
1999|), the surveys discussed above are not inconsis- 



tent out our principal results, which are: (1) the fraction 
of K-l-A galaxies (as defined in § [3.2^ in the field and X- 
ray luminous clusters at z « 0.3 is less than 5%; (2) it 
has not been shown conclusively that the z w 0.3 K-l-A 
fraction is significantly larger than the fraction at z « 0.1; 
and (3) there is no strong evidence yet that K-l-A galax- 
ies are more common in cluster environments than in the 
field. We suggest that the higher fractions of K-fA galax- 
ies found by the MORPHS group may come from the very 
different methodologies used, especially in the selection of 
the galaxy sample. We suggest that the simpler CNOCl 
sample and weighting functions are likely to produce a fair 
sampling of the cluster populations. In addition, there 
may be an actual difference in the populations of our two 
differently selected cluster samples. In that case, it re- 
mains to be seen whether clusters at these redshifts are 
better typified by the X-ray luminous CNOCl sample, or 
the more heterogeneous MORPHS sample. 

6. DISCUSSION 



iSpQj. reference, 7.6itl.l% of cluster galaxies, 
Wo{OII)< 5A (uncorrected for scatter). 



It was shown in Balogh et al. (1997) that the star for- 
mation rate (determined from [O II] emission lines) at all 

and 9.1itl.6% of field galaxies in our luminosity limit sample have Wo{HS)> 3A, and 
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radii within CNOCl clusters is always less than that in 
the field galaxy population. However, if a substantial frac- 
tion of cluster galaxies have undergone massive but very 
short starbursts, as postulated for example by Barger et 
al. (1996), the short duration of these bursts (salOO Myr) 
would mean that only a small fraction of cluster galaxies 
would be observed in the burst phase at a given epoch. 

According to the results of the PEGASE models, follow- 
ing a short burst a galaxy will appear as a FSB or PSF type 
for between 0.75 and 1 Gyr. Thus, the starburst rate can 
be determined from the abundance of these "remnants" . 
For our maximal sample (§2.3.3), the fraction of these 
galaxies in the cluster is less than 6.2±G.7%; this is an 
upper limit because it does not include the correction for 
non-uniform scatter. This implies that, in the last 2 Gyr, 
less than approximately 15±2% of the galaxy population 
may have undergone a starburst. Using our best estimate 
for the scatter correction from Table ^, the corrected frac- 
tion of cluster FSB and PSF galaxies is 2.2±0.9%, which 
reduces the estimated starburst population to only ~ 6% 
in the last 2 Gyr. This is much less than the 30% esti- 
mated by Barger et al. (1996). Part of this difference is 
due to the factor of two difference between the abundance 
of HDS galaxies in our respective samples, as shown in 
Figure E4 (significant at the ^ 2ct level). Secondly, our 
fractions are reduced due to the scatter correction, though 
Barger et al. include their uncertainties in their models, 
which should account for this effect if their error bars are 
reliable. Finally, Barger et al. actually model the HDS 
galax ies (which comprise 4.6±1.8% of our cluster sample. 



after iicatter correction) with post-starburst tracks; how- 
ever, some of these have [O IIJ emission (the A-|-em galax- 
ies), and cannot therefore be truly post-starburst. 

Poggianti et al. (1999) have recently suggested that the 
A-|-em galaxies are dusty starbursts which may later evolve 
into K-f A types; the strong Balmer lines arise because the 
OB-stars responsible for the emission lines are heavily ob- 
scured by dust, while the longer lived A stars have mi- 
grated out of their birthplaces and are more widespread 
throughout the galaxy, so that their light dominates the 
spectrum. In our cluster sample, we find 3.0±0.9% of the 
galaxies have A-|-em spectra, after scatter correction. This 
is almost two times larger than the K-t-A galaxy fraction 
and implies that the A-|-em phase must last at least twice 
as long as the K-l-A phase if the two types are evolution- 
arily linked, and if all A-|-em galaxies evolve into K-l-A 
types. Thus, these dusty starbursts should either be fairly 
long-lived, with lifetimes of more than 1.5 Gyr, or a large 
fraction of them never evolve into K-l-A types. The frac- 
tion of field galaxies classified A-|-em, 6.3 ± 2.1%, is twice 
as large as the cluster fraction. Thus, the cluster environ- 
ment does not appear to be responsible for preferentially 
generating these galaxy types, whether or not they are 
starbursts. It is also possible that, if this interpretation 
of A-|-em spectra is correct, some K-l-A galaxies are just 
a more extreme example, in which the [OIIJA3727 line is 
completely obscured, as pointed out by Poggianti et al. 
In this case, A-|-em and K-l-A galaxies are not evolution- 
ary counterparts, but representative of the same, starburst 
phenomenon. 

Ass uming that K+A galaxies are the result of recently 
termi: lated star formation, we can use the fractions of PSF 



and PSB galaxies in our sample to assess the relative con- 
tribution of starbursts and truncated star formation to this 
scenario. In general, PEGASE models of post-starburst 
galaxies spend roughly an equal amount of time in both 
PSB and PSF stages, while truncated star formation leads 
only to PSF galaxies, lasting about 300 Myr after the star 
formation activity ceases (see §3^). Thus, if starbursts 
are the dominant mechanism for generating HDS galax- 
ies, there should be (roughly) a n eq ual number of PSB 
and PSF galaxies. As shown in §4.3, the PSF fraction in 
both cluster and field environments may be much higher 
('^ 5 — 10 times) than the PSB fraction. If this overabun- 
dance is real, it suggests that at least some of the galaxies 
may have had their star formation truncated without a 
burst, as supported, for example, by recent simulations 
(Fujita et al. 1999). 

The lack of a marked excess of K-f-A or A-|-em galaxies 
relative to the field may suggest that the cluster environ- 
ment does not actively truncate or induce star formation, 
but merely inhibits it from regenerating once it has ceased. 
For example, in the galaxy formation model of Baugh, 
Cole & Frenk (1996), after a merger of two spiral galaxies 
succeeds in producing an elliptical galaxy, halo material 
may recollapse to form a new disk. This reformation of 
a stable disk could conceivably be inhibited in the clus- 
ter environment, due to the presence of strong tidal fields 
and galaxy harassment. Furthermore, the halo of gas sur- 
rounding a spiral galaxy, which may continually cool onto 
the disk, forming new stars, may be easily stripped away 



when a ga laxy falls into a cluster (Larson, Tinslcy & Cald- 



well 1980). The conditions for this to happen are much 



less stringent than those for stripping all the gas out of 
the galactic disk. In either of these cases, K-f A galaxies 
will generally evolve to become Passive within the clus- 
ter once they have exhausted their gas supply, whereas, 
in the field, many will again become star-forming galax- 
ies. Starbursts and K-l-A phases are then interpreted as 
natural stages in galaxy evolution (i.e., galaxies form stars 
in bursts, not continuously), and the cluster merely pre- 
vents a galaxy from recommencing its star formation once 
the activity has stopped of its own accord. We are cur- 
rently analyzing dark matter simulations of clusters to link 
a galaxy's star formation history to the time elapsed since 
it was accreted. Preliminary analysis shows that a gradual 
reduction of SFR following galaxy accretion adequately re- 
produces the radial gradient of star formation (Balogh et 
al., in preparation). 

There are several caveats which may alter the conclu- 
sions described above, and there are several ways in which 
future observations may address these issues; we list these 
caveats below: 

1. Once statistical correction weights are included, our 
spectroscopic sample is representative of a photomet- 
ric sample complete to Mr = —18.8-1-5 log h. When 
star formation is truncated in the more numerous, 
faint galaxy population, they will fade and many will 
drop below this limit. Indeed, the K-l-A galaxies in 
our sample appear to be of intrinsically low luminos- 
ity. It is also known that the less lumin ous galaxies 



evolve more significa ntly with redshift ( Lilly et al^ 
1995|; Lin et al. 1997), and thus contribute more to 



the Butcher-Ocmler effect. Thus, it is possible that 
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3here exists a large number of faint K+A galaxies 



which we are not detecting, and these may show an 
environmental preference. 

2. Whether or not an equal fraction of K+A galaxies 
in the cluster and field necessarily implies an equal 
production rate of starbursts in both environments 
is somewhat model-dependent, as it may only be a 
subset of the total population that is ever engaged in 
starburst activity. For example, the results of Ellis 
et al. (1997) and Barger et al. (1998) suggest that 
most cluster ellipticals formed from a single burst 
many Gyr ago, and have since evolved passively. If 
we assume that the Passive type galaxies play no 
part in cluster evolution, then the fraction of "ac- 
tive" galaxies in clusters that have a K+A spectrum 
is 6 ± 3%, (corrected for scatter effects). However, 
it is not yet clear how many of the Passive galax- 
ies are safely excluded in this manner; many of our 
Passive galaxies have spectra consistent with those 
of early type spirals and SO galaxies. Furthermore, 
it is not clear what fra ction of the field ga l axies are 
also "primo rdial" (e.g.. Kelson et al. 1997; Kodama| 
it al. 199"8|). 



et al. 1997). Work is currently being undertaken 



The [Oil] index is metallicity dependent, which af- 
fects star formation rates derived from its measure- 



ment; howeve r, this dependence is fairly weak (Ken- 



licutt 1992a) and unlikely to strongly affect our re- 
sults. More important is the fact that this blue fea- 
ture is quite sensitive to dust obscuration, and many 
of the PSB or PSF galaxies with Wo{OII)<^ SA may 
have significant, but dust-obscured, star formation 
(e.g., Poggianti et al. 1999). Infrared spectroscopy 
and imaging, and measurements of Ha emission, will 
be helpful in constraining this effect. 

The K+A popu lation is known to be morphologically 
heterogeneous ( Couch et al. 1998 ; Caldwell, Rose fc 
Dendy 1998 ) . If high resolution HST imaging reveals 
that K+A galaxies in clusters are morphologically 
distinct from those in the field, this may support 
models in which the mechanisms which are respon- 
sible for generating K+A galaxies are environment- 
dependent. 



5. Kauffmann (1995) and Andreon & Ettori (1999) 
have suggested that the X-ray luminous clusters ob- 
served at moderate redshifts are not the evolution- 
ary predecessors of the less luminous, low redshift 
clusters to which they are often compared. In par- 
ticular, the CNOCl sample is composed of the most 
X-ray luminous clusters in the universe and, hence, 
may have unique or unusual galaxy populations; this 
may be responsible for some of the discrepancy be- 
tween our results and those of Dressier et al. (1999). 
We note, however, that Andreon & Ettori find no 
dependence of the blue galaxy fraction on X-ray lu- 
minosity, specifically. 

6. Finally, the field sample drawn from the foreground 
and background of these rich clusters may be un- 
usual in some way, though results are generally con- 
sistent with the field samples of other groups (On 



to address this issue using the CN0C2 field galaxy 
sample (EUingson et al., in preparation). 

7. SUMMARY AND CONCLUSIONS 

We have presented a detailed analysis of the spectral 
characteristics of a large sample of galaxies in X-ray lu- 
minous clusters between z=0.18 and z=0.55, and of an 
identically selected field galaxy sample. We focus on three 
spectral indices: D4000, which traces the old stellar pop- 
ulation; Wo{H6), which indicates the presence of A-type 
stars and is sensitive to star formation that took place up 
to 1 Gyr ago; and the Wo{OII) index, which indicates 
current star formation activity. We compare our data to 
the model predictions of Fioc & Rocca-Volmerange (1997) 
and Bruzual & Chariot (1993), and draw the following con- 
clusions: 

• The radial trends within the cluster sample are con- 
sistent with a continuous age sequence, in the sense 
that the last episode of star formation occurred more 
recently for the outermost galaxies than for the cen- 
tral galaxies. 

• We define K+A galaxies as those with Wo{OII)< 
bA, and Wo{HS)> bA. Our measured K+A galaxy 
fraction is 4.4 ± 0.7% in the cluster and 1.9 ± 0.6% 
in the field, but this is only an upper limit due to 
the large index errors which tend to overestimate 
these numbers. Attempting to correct for this ef- 
fect, we find the true fraction of K+A galaxies to 
be only 2.1+0.7% in the cluster environment, and 
0.1+0.7% in the field. For our luminosity limited 
sample (galaxies brighter than Mr = —18.8+5 log /i), 
these corrected values are 1.5+0.8% (cluster) and 
1.2+0.8% (field). Our field results are consistent 
with the LC RS z = 0.1 fraction of 0.30% ( |Zabludofi 
et al. 199^ ) at the -- la level. 



• The fraction of cluster galaxies which are undergo- 
ing, or have recently undergone, a short burst of star 
formation is not significantly greater than the field 
fraction, at any redshift or distance from the cluster 
center. From the fraction of PSB and PSF galaxies, 
we conclude that less than '^10% of the galaxies in 
both the cluster and the field may have undergone 
short starbursts in the last 2 Gyr. 

• If all the A+em galaxies are dusty starbursts and 
progenitors of the K+A galaxies, then these must be 
undergoing fairly long-lived episodes of star forma- 
tion (> 1.5 Gyr). These galaxies are two times more 
common in our field sample than in the cluster. 

• We find that PSF galaxies may outnumber PSB 
galaxies by a factor of 5-10, in both cluster and field 
environments (though the samples are small) . Com- 
parison with models suggests that truncated star for- 
mation without a short starburst phase may play a 
significant role in galaxy evolution. 

• More photometry of the CNOCl fields is still needed; 
particularly K-band to constrain the total stellar 
mass, and U-band to independently measure star 
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formation rates. Infrared imaging and spectroscopy 
is needed to determine the amount of dust-obscured 
star formation, and HST images are required to sub- 
classify the various galaxies types morphologically. 

Keeping in mind the caveats listed at the end of §o, we can 
discuss some of the implications of our results. The gen- 
eral trend for higher redshift clusters to be bluer (the B-0 
effect) has not yet been linked to a population of galaxies 
unique to the cluster environment. The simplest expla- 
nation consistent with the current data is that the B-0 
effect in clusters largely reflects the increased level of star 
formation in the field at larger redshifts. Secondly, there 
is no evidence that recent (<1 Gyr) cluster-induced star 
formation is responsible for driving the differential evo- 
lution between cluster and field within the virial radius, 
since there is no tell-tale population of starburst or post- 
starburst cluster galaxies that does not exist in the field. 
We speculate that the presence of bursting galaxies, and 



the possible increase of such galaxies with redshift, reflects 
the general nature of star-forming galaxies, independent 
of environment. However, galaxies in dense environments 
are eventually prevented from bursting, perhaps because 
their halo gas reservoir has been stripped. This, coupled 
with a population of "primordial" ellipticals and perhaps 
some truncation due to ram pressure stripping of disk gas, 
could be responsible for the older stellar populations which 
inhabit cluster environments. 

The data used in this paper form part of the CNOCl 
study of intermediate-redshift clusters. We are grateful 
to all the consortium members and to the CFHT staff for 
their contributions to this project. MLB is supported by 
a Natural Sciences and Engineering Research Council of 
Canada (NSERC) research grant to C. J. Pritchet and an 
NSERC postgraduate scholarship. MLB is grateful to J. 
F. Navarro and A. Zabludoff for useful discussions. 



APPENDIX 



A DESCRIPTION OF THE PEGASE MODELS 



In this Appendix we will discuss in more detail the PEGASE models shown in Figure O, on which our galaxy classifi- 
cation system (§3.3) is partly based. 
The large, open star in Figure 



11 



represents the model result for a galaxy with a constant star formation rate, of any 
intensity. It represents the maximum US absorption one can expe ct in a normal galaxy; emission filling will reduce its 
value by between roughly 1 and 3. 5 A (Barbaro & Poggianti 1997). This point lies blueward of our D4000 cut, within 
the SF region. All of the galaxies in the SF region are either spirals, irregulars, or "starburst" galaxies, and d o n ot have 
Wo{H6) stronger than the model constant star-formation point (with the exception of NGC 3034, discussed in §3.2). Only 
if the star formation activity has begun recently (within the last 200 Myr, Barbaro & Poggianti, 1997) will US emission 
overwhelm the absorption and t he g alaxies will then occupy the short starburst (SSB) region of this figure. 

The long-dashed line in Figure O traces the evolutionary path of a galaxy with an exponentially decaying star formation 
rate; it evolves from "birth" at D4000« 1, redward. In this case, a decay time of r = 2 Gyr is adopted, and the galaxy 
is evolved to 11 Gyr, which we expect to be the maximum age of CNOCl galaxies, in reasonable cosmological modelsR 
Most of the local spirals in the Kennicutt (1992b) and Kinney et al. (1996) samples lie close to this curve. Galaxies of 
this type with ages greater than about 5 Gyr will lie in the Passive region; clearly this includes not only E and SO type 
galaxies, but also early type spirals. Note that the oldest models have Wo{HS)< 0; this reflects the nature of our specific 
index definition, which in this case is telling more about features in the continuum, rather than the intrinsic US absorption 
(see §m). 

The evolution of a galaxy produced by an initial burst of star formation lasting 200 Myr is shown as the short-dashed 
line in Figure O. After the burst ends, this spectrum always has stronger Wo{HS) than the exponentially decaying 
star-formation model, until D4000> 1.7. It spends roughly 300 Myr in the bHDS region and another 300 Myr in the 
rHDS region. 

The dotted line represents the evolution of a galaxy which underwent constant star formation for 4 Gyr, followed by a 
200 Myr starburst involving 30% of its mass, after which all star formation was terminated. This model track only follows 
the evolution after the end of the starburst. The model closely traces the initial burst model, but does not reach as strong 
values of Wo (HS) . The duration and strength of the burst are quite arbitrary; longer bursts and bursts involving smaller 
amounts of stellar mass result in weaker maximum Wo{HS) indices; this parameter space has been explored extensively 
by Poggianti & Barbaro (1996), though we warn that our Wo{H5) indices are defined differently and are not directly 
comparable with theirs. 

The rHDS galaxies can originate not only from a terminated SSB phase, but also from a galaxy in which long-term 
star formation has recently ended (e.g., pouch fc Sharpies 1987 , Newberry, Boroson & Kirshner 1990). We model this 
latter case by terminating star formation in a galaxy which has undergone constant star formation for 4 Gyr; this is 
shown in Figure n^ as the solid line. Once the star formation is terminated, the galaxy quickly reddens, passing into the 
rHDS region within 100 Myr. It then spends about 300 Myr in the rHDS region before Wo{HS) further weakens and the 
spectrum becomes a Passive type. 



^Beyond an age of 11 Gyr, the models only increase in D4000 at a rate of about 0.05/Gyr. 
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